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Abstract 

This paper studies a number of matrix models of size n and the associated Markov 
chains for the eigenvalues of the models for consecutive n's. They are consecutive principal 
minors for two of the models, GUE with external source and the multiple Laguerre matrix 
model, and merely properly defined consecutive matrices for the third one, the Jacobi-Pifieiro 
model; nevertheless the eigenvalues of the consecutive models all interlace. We show: (i) 
For each of those finite models, we give the transition probability of the associated Markov 
chain and the joint distribution of the entire interlacing set of eigenvalues; we show this is 
a determinantal point process whose extended kernels share many common features, (ii) 
To each of these models and their set of eigenvalues, we associate a last-passage percolation 
model, either finite percolation or percolation along an infinite strip of finite width, yielding 
a precise relationship between the last passage times and the eigenvalues, (iii) Finally it 
is shown that for appropriate choices of exponential distribution on the percolation, with 
very small means, the rescaled last passage times lead to the Pcarcey process; this should 
connect the Pearcey statistics with random directed polymers. 

1 Introduction 

In this paper we are concerned with the (generalized) minor processes associated to random 
matrix models that are related to very classical orthogonal polynomials. We also show their 
relation to last-passage percolation models. This project was partially motivated by the follow- 
ing open problem: finding a continuum random directed polymer interpretation for the Pearcey 
process, and a corresponding stochastic heat equation, in the same way that the Airy process 
has a KPZ / stochastic heat equation / random polymer interpretation; see the work of Amir, 
Corwin and Quastel [3]. A first step in that direction is to show that the Pearcey process 
appears as a limit of a last-passage percolation model; this is done in the present work. 

The first minor process arising from random matrix is the Gaussian Unitary Ensemble 
(GUE) minor process defined and analyzed in [17| by Johansson and Nordenstam. Following this 
result, other minor processes of classical random matrices are obtained [12]. The minor process 
of Laguerre Unitary Ensemble, aka complex white Wishart ensemble, leads to the generalized 
Wishart ensemble that was conjectured in p] and solved in [TO]. See also [13], [H], [2] and [1] 
for other minor processes related to random matrices. 
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Analogous to complex Wishart ensemble (aka Laguerre Unitary Ensemble with external 
source) that is a generalization of the classical Laguerre Unitary Ensemble (LUE), the GUE 
with external source is a generalization to the classical Gaussian Unitary Ensemble. In this 
paper we consider the minor process associated to this matrix model. 

In the same spirit as the complex Wishart ensemble and GUE with external source, there 
is a new matrix model that generalizes the classical Jacobi Unitary Ensemble (JUE) that we 
denote as Jacobi-Piheiro ensemble, and another generalization of the LUE that we denote as 
multiple Laguerre ensemble (to be distinguished with the complex Wishart ensemble). 

We describe the minor processes in a systematic way as follows. In all cases below, Wn, Xn, Yn 
denote the matrix of the first n columns of the M x A^-matrix W, X (resp. the M' x A^-matrix 
Y) for 1 < n < A^, where M, M' , N are positive integers and we assume M > N and M' > A^. 
For square matrices Zq and A of size A^, (Zq + A)n denotes the n-th principal minor of the 
square matrix Zq + A, also for 1 < n < A^. For the definition of X, we need parameters Oi that 
are nonnegative integers satisfying 

1 + ai < 2 + Q2 < • • • < A^ + aAr. (1) 



The four (generalized) minor processes Consider the spectra A^^-*, . . . , A*-^) of the follow- 
ing consecutive matrices 5„, where A^") = (aS"\ . . . , A^"^) and aJ"^ are in increasing order. 

1. (GUE with external source) The n x n consecutive matrice^ (minors of Sn) 



Sn = (Zo + A)n with 



Zo = GUE(Af), 

A = diag(ai, . . .,aN) 



2. (Wishart) The n x n consecutive matrices (minors of Sn) 



5„ = W*Wn with ( 



W = M X N matrix, 

^W^,=Af{0,-^), (3) 



where oi, . . . , a^v are negative parameters. 
3. (Multiple Laguerre) The n x n consecutive matrices (minors of Sn) 

( X = M X N matrix, 
S. = X:X. wrth < ^^. =^^(0,1/2 



ax,, = AA(0, 1/2) ( iorl<^<J+a„ (4) 



and Xij = otherwise. 



4. (Jacobi-Piheiro) The n x n consecutive matrices 



XlXn 



X = M X N matrix, as before, 



Sn = y " " V With <; y = M' X A^ matrix, (5) 
A„A„ + r„ r„ I ^Y^^ ^ ^^^^ ^^2), Qy.^. = aA(0, 1/2). 

Remark 1. Note that in the definitions of the first three minor processes, Sn-i is a minor of Sn, 
so the name "minor process" is assigned. Although in the last process, Sn-i is not a minor of 
Sn, by definition the numerator and denominator of 5„_i are minors of those of Sn respectively, 
and also in that case the eigenvalues of Sn-i and those of 5„ are interlaced as for eigenvalues 
of a Hermitian matrix and its minor. Thus the name minor process is also justified. 



^GUE(A^) is standardly defined as Hermitian matrices with independent normal entries, with 5R(Zo)ij 
AA(0,i), '^{Zo^j =AA(0, I) and {Zo).. ^ M{Q,1). 
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Remark 2. The Whishart minor process is a special case of the so called generalized Wishart 
random-matrix process |10j . and its properties has been already known. We include it in this 
paper for completeness. 

Previously studies minor processes are shown to be the continuous limits of special Schur 
processes and are equivalent to continuous last-passage percolation models with properly chosen 
parameters and possibly taking limit (see e.g. [9], [15], [10] and [12] for the cases most close to 
ours). The minor processes considered in our paper also have this property. 



(m, 1) 



(1,1) 



;i,^) 



Figure 1: Two up-right paths in n(m, n) that do not intersect. 



Consider the percolation model on the Z+ x Z"*" lattice {(i, j) | i,j = 1,2,...}, where on 
each site we associate a real weight Xij. Let n(m, n) be the set of up-right paths in the rectangle 
with vertices (1, 1), (l,n), (m, 1), (m, n), see Figure [ij Then we define the maximum of the total 
length of I non-intersecting up-right paths in that rectangle to be 



L^^\m, n) :-- 



E 



E 



max 

Pi,...,P;en(m,n) 
Pi,...,Pl non-intersecting'^"^ (hj)&Pk 



(6) 



In the case I = 1, L^^\m, n) is the length of the longest up-right path from (1, 1) to (m, n). Let 
pi,P2, ■ ■ ■ ;qi,q2, ■ ■ ■ be two sets of positive parameters, and let the Xjj's be independent and 
exponentially distributed of parameter VTjj = Pi + qj > 0, 



(xij >t) = e~ 



(7) 



Then the L^^\m, n) become random variables, which for fixed m, n are increasing as / increases. 
We define the n-variate random variables 



(m,n) 



(m,n) 



; • • • ) 



Where 5^/.^'")=lW( 



m, nl. 



(8) 



i=l 



(n) 

Now we state the distribution of the eigenvalues A^- in the minor processes defined above. 
Before the statement of the theorem, we first define the notation /u ^ of interlacing between 
an (n — l)-variate random variable fj, = {fii, . . . and an n-variate random variable i' = 

(z^i, . . . , f„) such that their components are in increasing order and are in a common domain /. 



fi ^ u if z/i < /ii < < /"2 < • • • < ^n-i < i^n, where fii, vi G /. 



(9) 
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Theorem 1. The interlacing sets of spectra X^^\ ■ ■ ■ , A^"), . . . , A*^^-* constitute an inhomogeneous 
Markov chain with transition probability: 



Pn-l,n(A("-l\A(")) 



1 



An(AW)nr=i^n(A! 



, .-(„-l) lA("-i)dAW' (10) 



where A„(A("')) = ni<j<A:<n('^i"^ ~ ^t^^^ Vandermonde determinant, the range I of the 

eigenvalues, the weight Wn{z), the function ip{z) and the constant Cn are given in Table [7| 



model 


/ 


Wn{z) 


Hz) 


Cn 


Pl(^) 


1. GUE ext source 

2. Wishart 

3. Mult. Laguerre 

4. Jacobi-Pineiro 


M 

[0,oo) 
[0,oo) 

[0,1] 


Z^n^-z 


1 

Z 

z 

z{l-z) 


2 

V27re 2 
(M-n)! 

a„! 

a„!(M'-n)! 
(q„+M')! 


■\/27V 

(M-1)!^ 
z"^ p-Z 
ail J 

{M'-iy. ^) ^ 



Table 1: The range /, the weight Wn{z)^ the function ip{z), the constant C„, and the probability 
density function pi {z) of the eigenvalue A^^^ of Si for the four minor processes. 

As a corollary to Theorem [l| we have the joint distribution of A^"^ in the minor processes. 

Corollary 1. In each of the four minor processes associated to the GUE with external source, 
Wishart ensemble, multiple Laguerre ensemble and Jacobi-Pineiro ensemble respectively, the 
joint distribution of the eigenvalues A^^-* G l'^ , G . . . , A^^^^ G is 



N 



N-l / n 



Wn{\ 



p(A«,...,A(-))=A.(A(-))n^-MA(-)) n n — ;:w;;:,wJ ^a(^ha(^... 

n=l \i=l Wn+liX] OV-lA- ') 



n=l 



or equivalently, (with the X^^i virtual variables introduced for convenience) 

N ^ N-l 



P(AW, . . . , aW) = A^(aW) n 7^-iv(A(^)) n det{Mxt^\x^^'^))^^l,, 
where </'n(Aj-"'\ A^"^^^) is given by 



(11) 



(12) 



n=l 



n=l 



and otherwise 



Wn+l[,x)'lp[x) 



<An(A^i,x) = l. 



(13) 
(14) 



Here the domain I, constant Cn and functions Wn,ip are defined in Table^ 

To prove Corollary 1 we use the Markovian property of X^"\ apply inductively, and 

(1) — 

note the initial condition that the probability density function of X\ in the four minor processes 



IS given bypi(z) in Table [!} O ne can check that the determinant in (|12l) encodes the interlacing 
property. 

The next theorem shows the equivalence between the minor processes and the last-passage 
percolation models. 



4 



Theorem 2. To the exponentially distributed percolation model n^.A^ with parameter 

TTij ■= Pi + Qj forl<i<M and 1 <j <N, with N < M, (15) 

we associate the variables /u -^'"^ , 1 < i < n as in 1^. Then the following holds: 

1. ForTTij = l—aj/\/^, the scaling limit . . . , u'^^^ of the percolation variables fi^^^'^\ . . . , 

vr' := lim ^ — , l<i<n< N, (16) 

M^oo Jm - - - , 



has the same joint distribution as X^^\ . . . , A^^-* in the GUE with external source minor 
process defined in 

2. (A special case of lT^ Theorem l.lj) ForTTij = —aj, the percolation variables iJ,^^^'^\ . . . , ^(*^'^) 
have the same joint distribution as X^^\ ■ ■ ■ , A^^-* in the Wishart minor process defined in 

3. ForTTij = i+ctj, the scaling limit . . . , v^^^ of the percolation variables . . . , 

v\ ' := lim Me-f"^ , l<i<n<N, (17) 

Af— 5>oo 

has the same joint distribution as the X^^\ ■ ■ ■ , \^^^ in the multiple Laguerre minor process 
defined in (Q. 

4- For TTij = i + aj, the variables v^^^ . . . , v^^^ obtained by exponentiating the percolation 
variables ^(^^''i), . . . , ^(*^''^), 



in) ._ 



e 



(M',n) 



l<i<n<N<M\ (18) 



have the same joint distribution as the \^'^\ ■ ■ ■ , A^^-* in the Jacobi-PiHeiro minor process 
defined in 

(n) 

Furthermore, we show that the correlation function of A^ in each minor process has a 
determinantal formula, and the correlation kernel has a double contour integral formula. 

Theorem 3. In each of the four minor processes associated to the GUE with external source, 
Wishart ensemble, multiple Laguerre ensemble and J acobi- Pineiro ensemble respectively, the 
correlation kernel of eigenvalues in Sm and Sn2 is given by 

(a) GUE with external source: 



-1 r Qiy-x)-w 



K{ni,x;n2,y) = 7^ ^ ™ 7;;^ — —dwl^^yln^^ni 



where the contour Fq encloses all poles of the form in the intergrand, and S = C + iM t 
lies to the right ofVa- 
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(b) Wishart (a special case o/ |3, Formula (15)]): 



K{ni,x;x2,y) 



n ■ T TT"2 t N ""^■^a;<j/-'i-ni<n2 



+ 



{2m) 



dz i dwe-'^^+y'" 



w — ai) w — z 



(20) 



where the contour Ta encloses all poles of the form Oj in the intergrand, and S encloses 
0, and does not cross or contain Ta- 



(c) Multiple Laguerre: 



K{ni,x;n2,y) 



-1 



o • r Tr«2 ? rdwl.x<y^ni<n2 



+ 



1 



dz (h dw 



X 



-2— 1„,UI 



{2mY Jr. '^"^ (.z - w)T{w + 1) nrii(^ - «^ 



(21) 



where Ta is a contour enclosing ai, . . . ,0^2; o,nd T, is a deformed Hankel contour going 
counterclockwise from —00 to —00 that encloses poles z = —1, —2, . . . (see Figure^ and 
the contour Ta ■ 

(d) J acobi- Pineiro: 
K{ni,x;n2,y) 



-dwl.x<y^ni<n2 



+ 



{2TTif 



dz i) dw 



X 



-z—l„,w 



T{w + M' + l)T{z + 1) Y[ll^{z - at] 
w T{z + M' + l)T{w + l)\Ylli{w-ai) 



(22) 



where the contour Ta is a contour inclosing ai, . . . ,an2, o,nd T, is a contour going coun- 
terclockwise enclosing — 1, — 2, . . . , —M' and the contour Ta- 




Figure 2: The deformed Hankel contour. It comes from —00, ends at —00 and enclose (— oo,0) 
counterclockwise. It keeps a constant distance to the real line at —00. 



Remark 3. As a special case of the GUE with external source minor process, if all = 0, our 
model becomes the well known GUE minor model. In |lll Formula (21)], the correlation kernel 
of the GUE minor model was obtained as the sum of two contour integrals. Due to a different 



choice of scaling convention, the kernel in |1H Formula (21)] is related to our kernel in (19) by 



i^F-F''(6,m; 6,^2)^6 = 2^'''-'''^/^V2Kinux;n2,y)dy 



x=V2i,y=V2rj 



(23) 
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The correlation kernel of the GUE minor model was first discovered in |17l Definition 1.2], (see 
|18l Definition 1.1] for erratum) by Johansson and Nordenstam in terms of Hermite polynomials. 
Their kernel is related to ours by 



K3:^^ir,(;s,ri)d£, = e 4° K{ni,x;n2,y)dy 



(24) 

ni=s,n2=r ^ ' 



Although the formula of the correlation kernel in our paper is slightly different from those in 
previous literature, they define the same correlation function of eigenvalues of minors, since the 
difference is simply a conjugation and change of variables. Hence Theorem [T] is a generalization 
of previous results. 

Remark 4. If we let all On = in the GUE with external source minor process or all an = M — n 
in the multiple Laguerre minor process and the Jacobi-Piheiro minor process, they are reduced 
to the GUE, LUE and JUE minor processes respectively, and we can check that the correlation 
functions in the special cases agree with those obtained in [12j for these three minor processes. 

It is not a coincidence that these four minor processes share so many similarities. They 
are closely related to the so called very classical multiple orthogonal polynomials, namely the 
multiple Hermite (to GUE with external source), multiple Laguerre of the second kind (to 
Wishart), multiple Laguerre of the first kind (to multiple Laguerre) and Jacobi-Piheiro (to 
Jacobi-Piheiro) respectively. Here we state without proof the correlation kernel in the Wishart 
minor process (|20|) can be written in the form of multiple Laguerre polynomials of the second 



kind when ni > n2 (see (234) and (235) for the definition of -P(a„j,a„j_i,...,a„^^j;M-rii)(3;) and 
Q(a„2 ,ari2-i,---,a.nf. ■,M—n2 ){y)) 



K{ni,x;n2,y) — (-1)"^ -^(Qui ,a„-, -i,...,an,„ , -, ;M-ni)(^)Q(an,,a„,_i,...,ant.;M-n2)(?/)- 



(25) 

When m < n2, the kernel is also related to multiple Laguerre polynomials of the second kind, 



but the relation is not so simple. Although we do not prove (25), we indicate that it can be 
proven based on the joint probability density (12) and the algebraic result in Lemma [2] of the 
next section, leading to (20), and we prove similar formulas for the other three minor processes 
as intermediate steps in the derivation of double contour integral formulas (19), ( |21[ ) and ( |22[ ). 

It is well known that in the Gaussian unitary ensemble, Laguerre unitary ensemble and 
Jacobi unitary ensemble, if the dimension approaches infinity and we consider the correlation 
among eigenvalues at the edge of the support of their limiting empirical distribution, then we find 
the correlation kernel has the limit as the Airy kernel. If we consider the corresponding minor 
processes, it is shown in [12] that the joint distributions of eigenvalues of consecutive minors 
around the edge of the limiting empirical distribution has the limit as the extended Airy kernel 
that defines the Airy process. Since the minor processes discussed here are generalizations to the 
three minor processes analyzed in |12j . one may expect that they can realize more complicated 
correlation kernels as their limiting kernels as the dimension of minors approaches infinity and 
the parameters Oj or are chosen properly. Indeed, in |8], Borodin and Peche shows that the 
generalized Airy kernel with two sets of parameters can be realized in the limit of the generalized 
Wishart random-matrix process, which is a generalization of the Wishart minor process in our 
paper. In this paper, we show that if we choose parameters properly, the Pearcey kernel can be 
realized as the limit of the correlation kernel of the multiple Laguerre minor process. Similar 
result should hold for the other three minor processes considered in this paper, but we only 
analyze the multiple Laguerre case for brevity. 
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We consider one special case of the multiple Laguerre ensemble, which is the analytically 
most feasible one (besides the white Wishart ensemble), such that the parameter ai are of only 
two values, namely, half of them are na and the other half nb. When Oj are chosen in this 
simplest way, however, the region that the nonzero entries in X occupy is a composition of 
two trapezoids, aesthetically not of the simplest shape, compared with the composition of two 
rectangles. In the later case we also observe the Pearcey process when taking limit properly, 
but we omit the details for brevity. 



The trapezoidal multiple Laguerre minor process We define the M x N random matrix 
X as follows, depending on a large integer n and two parameters a, b. Let the left n columns of 
X be determined by the parameter a and the other N — n columns by 6, such that in the i-th 
column where i < n, the top i + [na\ entries are in i.i.d. standard complex normal distribution 
and the bottom M — i— [na\ entries are zero, while if z > n, the top i + [nb\ entries are in i.i.d. 
standard normal distribution and the bottom M — i — [nb\ entries are zero. We shall take the 
limit n — )• oo, and assume that is large enough, (say, at least greater than 2n), and M > N 
is large enough so that M >i + na and M >i + nb va. any column. In terms of parameters ai , 
the M X N random matrix X is characterized by 



Qi = a2 



Oir. 



\na\ , and a. 



n+l 



an+2 



[nb\. 



(26) 



For our purpose to analyze the limiting Pearcey process, a and b are chosen in the way that 
(see Figure [s]) 



the graphs of f{x) 



+ 



{x — a)^ (x — b) 



and g{x) 



1 

X 



intersect at a unique point xq on the interval (a, b). 



(27) 







Figure 3: An example of a, b and xq, where 
a ^ 1.8748, b ^ 8.0752 and xq ^ 4.6305. 



Figure 4: The central part of F^, S°^,rg^ 
and FajSjFf, around respectively, where 
poo ^ poo u poo (^ggp_ f = fa U ffe). For 
poo ^ Y,oo ^ poo ^ ^Yie rays go to infinity, while 
for ra,S,rf,, the contours extend to the 
ends having magnitude n*. 



Theorem 4. Suppose a is large enough and b is determinted by a by (27). Upon the change of 
scaling 

ni = [n(2 + cin^ 2 s)J , ^2 = [n{2 + cin^a t)J , 



x = cue, 



(1 — c^n iu), y = cnc2 [1 — c^n 'iv), 



(28) 
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where the constant c is expressed as 



.1 In 
c = rEoexp( \ -), 

xq — a xq — 



(29) 



and the constants ci,C2,C3 depending on a,b,XQ are defined in (177), the correlation kernel 



(21) of the multiple Laguerre minor process with parameters ai defined hy a,b via (26), upon 

s, 

K{ni,x; 712, y)dy = K^{s, u; t, v)dv, 



conjugation, becomes the extended Pearcey kernel, i.e., for fixed s,x,t,y 

x"^o((a;o-6)n)-"i 



lim , 

n-5-oo 2/"^o((xo — ojnj~"2 



(30) 



with 



K^{s, u; t, v) 



1 



dz 



dw 



e 4 



1 



1 



t>s 



(27rz)^ Jfoo e"x+^~"^ w — z 



^/Mt 



_g 2{t-s) 



(31) 



where the contours T°° = U and are shown in Figure U 



The condition that a is large enough is technical. See the discussion in Remark 11 



Outline of the paper 

In this paper, the joint distribution function of the eigenvalues in the minor processes are derived 



in Sections 2.1 3.1 and 4.1 using the idea of corank 1 projection used by Forrester et al. in [12] 



and |14j . The derivation of the determinantal kernel in Sections 2.2, 3.2 and 4.2 from the joint 
distribution function is based on Lemma [2] that was obtained by Borodin, Ferrari, Prahofer 
and Sasamoto in [71 Lemma 3.4]. Then we do the asymptotic analysis for a special case of the 
multiple Laguerre minor process to show the occurrence of the Pearcey process as the limit. 
The relation between the minor processes and percolation models and Schur processes via RSK 
correspondence in Section [g] follows the argument by Forrester et al. in similar models, see |12j 
and |14l Appendix A]. In Appendix [A] we derive new contour integral formulas for the two kinds 
of very classical multiple orthogonal polynomials, namely the multiple Laguerre polynomials of 
the first kind and the Jacobi-Piheiro polynomials, which do not appear in literature according 
to our limited knowledge. The construction of random matrix models related to the two very 
classical multiple polynomials brings new types of random matrix models, and the Pearcey 
kernel is seen in the minor processes for the first time. 
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2 Joint distribution of eigenvalues of minors in GUE with ex- 
ternal source 

2.1 Transition probability of the Markov chain A*^"^ 

The results in this subsection depend on the following technical lemma. 

Lemma 1. Suppose H' is an (n — 1) x — 1) fixed Hermitian matrix with distinct eigenvalues 
^ = (fii, . . . , fJ-n-i) in increasing order. Let H be the n x n random Hermitian matrix defined 
by 
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• The upper-left (n — 1) minor of H is equal to H' . 

• Denote the {n — 1) dimensional column vector h := ■ ■ ■ , ^n-i.n)"*"; the last column 
of H without the last component. The components of h are independent complex random 
variables in standard normal distribution, i.e., p{^hin) = Af{0, ^) and p{^hin) = A^(0, 

• hnn, the lower-right entry of H, is a real random variable independent of kin (i <n), and 
it is in normal distribution: p{hnn) = A/'(a, 1), where a is a real constant. 

Then the distribution of eigenvalues of H , denoted by X = (Ai,...,A„) with Ai < ••• < A^, 
satisfies the interlacing property /i ^ A given by ^ and their joint distribution is 



p^(A) = CA„(A)e^'^'=^ ^ "l^^A, where C = -— . (32) 

V27r A„_i(/i) 



Remark 5. The distribution of eigenvalues of H depends only on the eigenvalues of H' . 

Proof. The proof of Lemma [T] is analogous to the proof of ^ Section 4]. Suppose U G U(n — 1) 
is an (n — 1) dimensional unitary matrix depending on H' , such that UH'U~^ = diag(//). Then 
the n X n random Hermitian matrix 

H = (^^'""ll^^ ) , where h := (/ii, . . . , hn-i) = Uh, (33) 

is conjugate to H, and they have the same eigenvalues and the same characteristic polynomials 
Ph{z) = pfj{z). Note that components of h are also independent random variables in standard 
complex normal distribution. By direct computation, we find the relations between the charac- 
teri 
H': 



teristic polynomials ph{z) = Pg{z) = YVi=i{z — Aj) of H and H and Ph'{z) = Y\A=ii^ ~ f^i) 



Ph{z) hi^ , . 
]-::=z-hnn-y , 34) 

PH'{z) ~{^-l^i 

= , _ (^^(A) - - -{al{\)a^{^,) + ct2(/.) - ct2(A) - a?(/z)) + 0{\), (35) 

PH'{.Z) Z Z^ 



where cti(A) = o"2(A) = ^i^jXiXj are elementary symmetric polynomials. If \hi\ are 

all positive, the roots of ph{X), which are eigenvalues of H and H, satisfy the strict interlacing 
condition X ~< ^, which means 

Ai < ^1 < A2 < • • • < /in-l < An- (36) 

To see that, we notice the limiting behavior oi ph{z)/ph'{z) 

r PhW , y Ph{X) , . 

hm = ±00, hm = =Foo, (37) 

A^±«D ph' (A) x^(^l^)± Ph' (A) 

and find that there is one root in each of the n intervals with ^j, ±00 as endpoints. 

Given \hj\'^ G M+ and hnn G there is a unique A = (Ai, . . . , A„) satisfying (36). On the 
other hand, given A satisfying (36), there is a unique array of G M^,/i„„ G M such that 
Ai, . . . , A„ are roots of ph{z). To see that, we identify the residues at z = fij of the right-hand 
side of (34) with that of the left-hand side, and obtain 

,7 ,2 PnilJ-i) . . . .„„^ 

= T^^' ^ = !,■••, 38) 
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then we identify the constant terms in the right-hand sides of (34) and (35), and obtain 



Kn = 0-1 (A) - cri{fi). 



(39) 



The argument above also shows that the probabihty that some eigenvalues of H coincide with 
eigenvalues of H' is the same as the probability that some hi = 0, which is 0. Thus in order to 

^.1 



find the probability density of A, we need only to find the probability density of p G M+ and 



hnn G and the Jacobian determinant of the map from A to /inn, since 



p(A) = . . . , |/l„-l|^, Kn) 



d{\hi\'^, . . . , |/ln-lP, hnn) 



d{X) 



By (38) and (39) we obtain 



di\h\\..., 



d\hj\' 



hn—l\ 1 hnn) 



Aj — fj-j 



n-l 



dhr 



dX, 



d{X) 



Yl\hj\'^ det{C) = (-1) 



-1 An(A) 

An-l(/W)' 



where 



C 



/ 1 



V 1 



1 \ 



1 / 



Here we use the formula (38) of \hi\ and 



det(C) = (-l)^=-^ An(A)An_l(^) 



nr=in"=i'(A.-/i.) 



from the Cauchy determinant formula. 

On the other hand, by the distributions of and h^, 



(40) 



(41) 
(42) 



(43) 



(44) 



n— 1 1 5 hnn) 



2tt 



27r 



e-l EILi A?+| Eir/ M?+a(EILi a.-EIT/ m,)- V 



(45) 



where the second identity follows from identifying the —'^^=i\hi\'^ with the residue of the 
right-hand side of (34) at oo, and calculating the residue from the right-hand side of (35). Thus 
Lemma [T] is proved. □ 

Proof of Theorem^ ( GUE with external source). To prove that the random variables A^"'^^ . . . , A*^") 
constitute the Markov chain with transition density given in (10), it suffices to show the identity 
of limiting conditional probability density of A*-"-* that 



lim P(A(") I A« G /1(/^^'^),---,A("-i) G /r'(/x^""'^)) = ^„(.-i)(A(")), 



(46) 



where P.,(n-i) (A(")) is defined in ([32]) and 



/J(^{J)) = {{xi,...,Xj) \xi<--- < Xj and \xi - //^l < e}. 



(47) 
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Without loss of generality, here we assume that for all j = 1, . . . , n — 1, fi^ , . . . , ji - are distinct 

real numbers in increasing order. Under the condition A^-'^ G l| (/i'--'-*), j = l,...,n — 1, we 
assume the conditional distribution of Sn-i, the n — 1 minor of the matrix Zq + A, is given by 
fe{Hn-i)dHn-i- The distribution function fe{Hn-i) satisfies a property that fe{Hn-i) = if 
the eigenvalues of Hn-i, ordered increasing ly, is not in /^-^(^u^"-!)). This is a direct consequence 
of A("-i) G /f-i(;u("-i)). 

Then using the general formula for conditional probability density functions 

fxix \Y£A) = J fx{x I y G A and Z = z)fz{z \ Y G A)dz, (48) 

we have 

P(AW I A(^^) G j = 1, . . . ,n - 1) = j P(A(") I 5„_i = Hn-i)fe{Hn-i)dHn-i. (49) 

Here P(A''"''* | 5„_i = Hn-i) denotes the conditional probability density function of A^"^ that 
are the eigenvalues of 5"^ (in increasing order), the n-minor of Zq + A, where the condition is 
that Sn-i is equal to Hn-i- The value of the density function P(A^"^ | = i^n-i) depends 
only on the values of the eigenvalues of Hn-i, as a consequence of Lemma [Tj From (32), we 



know that P{X^"'^ \ Sn-i = Hn-i) depends on the eigenvalues of Hn-i in a continuous way, 
given that the eigenvalues are distinct. Then we know that for any A^"-* and 5 > 0, there is an 
e such that if the eigenvalues of Hn-i is given by ^ G (//("- 1)), then 

|P(A(") I Sn-l = Hn-l) - P,(„-i)(A("))| < 5. (50) 



Then (49) and ( |50[ ) imply that if e is small enough for ( |50[ ) to hold, then 

|P(A(") I A(^-) G (;u(^-)), j = 1, . . . , n - 1) - P^(„-.)(A(-))| < S. (51) 

Taking the limit e — )• 0, we obtain (46). Thus we prove Theorem [l] in the GUE with external 
source case. □ 

2.2 Correlation kernel of the Markov chain A*^"^ 



For the ease of derivation in this subsection, we re-express the formula ( 12 ), the joint probability 
density of A^"^ : 



...,,.„,AS"\A5."+^)))^;i, det(*j^_,(Af )))f,, 

> \n=l / 

where the function 0„(A-"\ A^"""''^'') (i 7^ n + 1) has explicit formula 

</.„(x,2/) = e(""-'^"+i)"l,.<j„ (53) 

and 

^f{x) = e-^+"-^'P(,,,,,_,,...,,,_^^,)(a;) = -=: / e^— J](, - a^.,+i)(i., (54) 

and -P(aiv,ajv_i,...,a]v_,+i)(^) the i-th degree (monic) multiple Hermite polynomial of type II 
(see Appendix [a]) . 

In the proof of Theorem 5|[a) we first express the correlation kernel in terms of multiple 



Hermitian polynomials, and then write it in the form of double contour integral. 
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2.2.1 The correlation kernel and multiple Hermite polynomials 

The expression of the correlation kernel in terms of multiple Hermite polynomials, and even the 
fact that the Markov chain A^") is determinantal such that its properties are captured by the 
correlation kernel, is based on the following general lemma of Borodin, Ferrari, Prahofer and 
Sasamoto [7]. 

Below we state the lemma. First we fix / as a subset of M (in our paper, / is taken as 
(—00,00), [0,00) or [0, 1]). Suppose the functions a{x,y),b{x,y),c{x) are defined on / x / or /. 
Then the operator * used in the lemma is defined as 

a*b{x,y) := J a{x, z)b{z,y)dz, a*c{x) := J a{x, z)c{z)dz, c*a{x) = J c{z)a{z,x)dz. (55) 

Lemma 2 ([71 Lemma 3.4]). Suppose we have a signed measure on I x X • • • X given in 
the form 

N-l 

f- n det(</)„(A!"\ +^)))^+i, det(*^_,(Af )))f,=i, . • • , A^) E (56) 

n=l 

where (\>n{}\'^ ^ A^"^^^) stands for the value of the two-variable function (j)n{x, y) at x = A^-"'*, y = 
^("+1) iji<^ji otherwise the value of the one-variable function (pn{X^^i,x), and Zj\[ is a 
normalization constant. If / 0, then the correlation functions are determinantal. 
To write down the kernel we let 



^(ni,n2)(^ l*?^"! * (<^ni+l * i' ' ' * (Pn2-l))ix , y) Ui < 712, 

' I ni > 712, 



(57) 



and for 1 < n < N , 

^^,.(x):=0("'^)*^j^_^.(x), i = l,2,...,iV. (58) 
Set (j)Q{X^^\x) := 1. Suppose the functions defined on I 

</-o*</.(i'")(Af\x), (</>„_2*0("-i'"))(Ai"_-'\:E), ,/.„_i(A(r-'),x), (59) 

are well defined, then they are linearly independent and generate an n-dimensional space Vn- 
Define a set of functions <l>"(x) fj = 0, . . . , ri — Ij spanning Vn by the orthogonality relations 

= ^u, /orO<i,j <n-l. (60) 

Under Assumption A: 

<A„(ASi,a;) = c„<+')(x), 7i = l,...,iV_i (61) 
for some Cn 7^ 0, the kernel takes the simple form 

Kinux- 712, y) := -<^("^'"^) (x, y) + ^ Kl-ki^)Kl-kiy)- (62) 

k=l 

Remark 6. Our version of the lemma is slightly different from that in [71 Lemma 3.4], in the 
sense that they consider discrete random variables but we consider continuous ones supported 
on /. The proof in their paper can be used for Lemma [2] with little change. 
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We first compute the correlation kernel under the assumption that oi > 02 > • • • > cat by 
the application of Lemma[2]with / = (— cxd,oo). We need the following identity that if (f>n{x,y) 



is given by (53), (pni^l^^i,^) = 1 and "^%_j{x) is given by (54), then 



dn-iX CLfioy 



:,iy-x)z 



2vri /r, n"ini+i(^--j; 
where Fa is a large enough contour enclosing all poles am+i, • • • , 0^2; 

andQ 



(63) 



^]{x) 



^dnX r 

\/27r« Jc 



27ri Jc+i 

,2 



6 2 \{s - an-k+i)ds 



27ri Jc+ 



(a„,a„_i,...,a„_j+i) l^ji 



e 2 



ds 



j = 0,1,. . . ,n - 1, 

i = -l,-2,...,n-A^, 



(64) 



where C is a large enough real number such that C + iM t lies to the right all possible poles of 
the integrand, and 



i>^i.-)=\ n 



1 



,fc=j+2 



aj+i — a/s 



o{a-i+l-o.n)x 



j = 0,l,...,n-2. 



(65) 



Note that from ([sol) and ([65|, the vector space T4 is spanned by e("i an)a::^g(a2 f^")^, . . . , e^"" '^"^^ 



and $^(x) are vectors in Vn defined by the orthonormality (60 ). Since ^'"(a;) are defined by mul- 
tiple Hermite polynomials of type II, $"(x) are hence defined by multiple Hermite polynomials 
of type I. Using Proposition [2] in Appendix \K\ we have 



$,"(x) = e^--"(^)Q(„„,„„_,,...,„„_^,)(x) 



27r27ri Jt 



+2 



dt 



ni=o(*-«^ 



(66) 



n—k ) 



where (5(a„,a„_i,...,a„_j)(^) is the multiple Hermite polynomial of type I. 



It is clear that the Assumption A is satisfied, as both functions (/)„(A^'^]^, x) and f^^^'^^' {x) in 



(61 



are constants. Therefore the correlation kernel is given by (62) with (j)^^'^''^^\x,y) expressed 

ni 



in (63), ^^^^_f^{x) expressed in (64) and ^^^^_f^{y) expressed in (66). 



Note that the correlation kernel K{ni, x; n2, y) is analytic in ai, . . . , cat, by (62 ). By analytic 
continuation we can remove the restriction ai > 02 > • • • > cat. 



Proof of (63), (64) and (65). We prove (63) by induction. If ni = 722 — 1, it is clear that (63) 



holds. If it holds for ni = n2 — k, k > 1, then for m = n2 — A; — 1, 



In the evaluation of ^q^\x) in (64 1 and below, we take the notational convention that n?=i(') = 1- 
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-oo 



,(y-w)z 

"2 — 7Z, — ^;^dzlLx<w<ydw 



-/r^ nn'+2(^-«i, 

yz-a„^+ix r-y 



„-^a;-a„2y r Jy-x)z _ ({/-a;)a„j_ 



27ri Jr„ nn?+i(^-a. 



27ri 



r, Wll+i{z-aj) 



27ri 7r, nn?+i(^-«i 



(67) 



where in the last step we use the vanishing of the contour integral jp 



„no —, sdz, which can 

n„i+i(^-aj) 

be seen by deforming Tq to a contour about oo, and using n2 — ni > 2. 

The proof of (64) is based on the contour integral formula (229) of multiple Hermite poly- 
nomials of type II in Appendix [A| From (58), (63), (54) and (229), we have, where C + iM t is 
to the right of Fq , that 



2i:i 



,{y-x)z 



-dz 



e 2 



+aNy 



e 2 ' (s — ai)dsdy 



N 



=j+i 



ds (h dz 



e 2 



n (s-ai) / e(^-^)^dy 

I 1 J X 



l=j+l 



-(27r)3/2 



N 



C+iRt 



dse'^ (s - «i) f ^^-2 



/27ri Jc+t: 



, sl_,x UZj+iis - ai) 
dse-2 — 



nfc=n+i(s-afc) 



(68) 



and from (68), (64) is easy to obtain. 

The proof of (65) is straightforward, at least if oi, . . . ,a„ are in strictly descending order. 
Using (t)n{^^n+i^ ^) ~ 1 ^^"^ (63), we conclude that 

^aj+iy-a„x f p{x-y)z 



'OO 

a„x 



27ri Jr, nfc=j+2(^-afc) 



dzdy 



dz 



27r« Jr. \[l=j+2{z - ak) 



^{a,+^-z)y^y 



^(ajj^i—an)x 



-1 



2^« Jr. nLi+i(2-afc) 

g(aj+i-a„)x 



(69) 



n 
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In the above derivation we assumed that the contour Ta encloses aj+2, 0^+3, . . . ,an but hes to 
the left of aj+i, which enables us to do the y integration and finally deform Ta to a large contour 
about 00, picking up the residue just at z = aj+i. The argument above relies on the assumption 



that aj+i > aj+2 > • • • > fln, but it is clear that (65) holds when the assumption is removed, 



due to the analytic continuation. 



□ 



2.2.2 Double contour integral formula of the correlation kernel 



From (68) and (66), 



"2 



ds (f> dte 2 



N 



(2vri)2 JC+. 
Using the identity that for any j < N 



2 -xs+yt I ifc=n2 + l ~ °fc) y> Ilk=j+1 ^k) 

UZm+lis - ai) Wi=j{i - a«) 



{s-t) 



n 



N 

k=j+l 



I (t - a-i) 



iN 



Y\lAt-ai) 



N 

n 

k=3 



ak 



and the telescoping trick, we find 
Uk=j+iis - ak) 



t - ak 



N 



N 

n 

k=j+i 



ak 



t — 



1 n£i(t-«i) 



n 



s — t-"--"- \t — ak 

k=l 



ak 



N 

n 



s — t \t — ak 

k=n2+l 



ak 



(70) 



(71) 



(72) 



Substituting (72) into (70), we obtain, letting the contour C + iM t of s be to the right of Fa, 

712 



EKl-ki-)Kl-kiy) 



k=l 



(27ri)2 



ds 



C+ii 



dte 2 



-xs+yt 



UkUi^-^k) 1 

U?Mt-ai) s-f 



(73) 



Hence we prove Theorem ^[a) from ( 63 ) , ( 73 ) and ( 62 ) in Lemma [2| and the fact that we may 
conjugate any Fredholm determinant without changing its value. 



3 Joint distribution of eigenvalues of minors in multiple LUE 

3.1 Transition probability of the Markov chain A*^"^ 

Analogous to Lemma [T] in Section [2j we need the following technical lemma in this subsection. 
In the statement of the lemma and later in this section, we denote /i, ^2, • • • be nondecreasing 
integers such that In > n, and denote the nonnegative integers a„ = — n. 

Lemma 3. Suppose X' is an In-i x ('^^ 1) fixed rectangular matrix with distinct singular values 
fj, = {y/JIi-, • • . , ^^n-i) where = (/ii, . . . , fin-i) is assumed to be < m < ■ ■ ■ < Hn-i- Let X 
be the In x n rectangular random matrix defined by 

• The upper-left In-i x (n — 1) block of X is equal to X' . 
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Denote the In dimensional column vector x := . . . , x^^^n)-*", the last column of X. 

Then components of x are independent complex random variables in standard normal 
distribution, i.e., = AA(0, 1/2) and p{Qxi^n) =AA(0, 1/2). 

All entries in the lower-left (In — In-i) x (n — 1) block of X are 0. 



Denote the singular values of X by \/Ai, . . . , v\i where A = (Ai,...,A„) is assumed to be 
< Ai < • • • < An. Then the joint distribution of A is 



n-l 



p^(A) = C7-iA„(A) n^^^"^' ^ = "n!A„_i(^) Yl //f +^e-'^\ (74) 



\i=l 



i=l 



Proof. For the convenience of the proof, we denote S = X*X and S' = (X')*X' , and note that 
the eigenvalues of 5 and 5' are Ai, . . . , A„ and //i, . . . , respectively. 

Performing the singular value decomposition, we find the unitary matrices U G U(/„_i) and 
V e U(n - 1) such that 



X' = U 












7*2 



\ 








V*. 



(75) 



V / 

Let {/ = [/ e I(i„-i„-i)x{i„-in-i) e U(/„) and y = y e /ixi G U(n), we then have 

X = uxv*, 

where 



(76) 



X 








V ••• 



xi,„ \ 




^n.n 







(77) 



and in the last column of X all components are in i.i.d. standard complex normal distribution. 
Below we use the fact that the eigenvalue distribution of S is the same of that of 



S := X*X 











\ 



■\/Jjm—lXn—l,n 



(78) 



Now we consider the relation between the characteristic polynomials ps{z) of S (that is the 



same as p^{z) of S) and ps'{z) of S", analogous to (34) and (35). By direct calculation, we find 
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that 



n— 1 



Psjz) ^ 

PS'{z) ^ ■^^-Mi 

M^) = , _ (^^(A) - - i(ai(A)ai(/i) + a2(^) - a2(A) - + 



where di, ct2 are defined the same as in (|35|), and in (|79|) 



i = 1, . . . , n — 1, 



I = n. 



(79) 
(80) 

(81) 



Noticing the hmiting behavior of P5(z)/p5"(z) analogous to (37) yields the strict interlacing 
condition A ^ /i when ^i, . . . ,Cn-i are all positive, that is, 



< Ai < /ii < A2 < • • • < ^„-i < A„. 



(82) 



On the other hand, given A satisfying ( |82[ ), there is a unique array of ^1, ... , ^n-i G ^ 
la 



To see that, by the calculations of residues like in (38) and (39), we have 
PsiiJ'i) 



]Xl=i{z-Xk) . . 
res % . z = l,...,n-l, 



= a,{\) - ai(^) -yii= res nfc=i(^ " + V^' ^es ^^'='\' ~ 



n-l 



(83) 
(84) 



i=l 



res 



Uk=li^-^k) _ 117=1 



zmziiz - ^^k) nr=iV/ 

Then we can express p{X) analogous to (40) 

d{^l,...,^n) 



5(A) 



Like (41), 



— ^ = - — — , for J = 1, . . . ,n - 1, and — 
dXi Xi- fij dXi 



n-l 



dXi 



(85) 



(86) 



and like (42), 



^n-l(/^) riLiVj' 



(87) 



where the n x n matrix C is defined in ( 43 ) , and in the last identity of ( 87 ) we use the formula 



On the other hand, by the definition (81 ) of ^1, . . . , they are independent, 2^„ is X2{a +1) 
in distribution and 2^i, . . . , 2^„_i are all xl in distribution, and so by expressions (83) and (84), 
(note On = In — n) 



p(6,...,en) = e-^"i«'^e-« 



1 /n-=iA, 



Substituting (87) and (88) in to (85), we prove Lemma 3 



□ 
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Proof of Theorem^ (multiple Laguerre case). Recall in the proof in the GUE with external 
source case, the Markov property of the process relies on the fact that the distribution of the 
eigenvalues of H in Lemma [T] depends on only the eigenvalues of H' , but not its eigenvectors. 
In Lemma [Sj we also have that the distribution of the eigenvalue of S" = X*X depends only on 
the eigenvalues of 5' = {X')*X', but not its eigenvectors. 

(7) (i) 

Like in the proof in the GUE with external source case, suppose for all j, , ■ ■ ■ ,^ are 



distinct real numbers in increasing order, and define as in ( |47| ). Then analogous to (49), 

P(A(") I A(^') G j = 1, . . . , n - 1) = j P(A(") I 5„_i = Hn^i)f,{Hn^,)dHn-u (89) 

where the distribution function f^[Hn-i) is defined in the same way as that in the proof in the 
GUE with external source case. Then after the same argument, we derive the limiting identity 
analogous to (46) 

P(A(") I A(^-) = j = 1, . . . ,n - 1) = P^(„-i)(A(")), (90) 
and we prove Theorem [T] in the multiple Laguerre case. □ 

3.2 Correlation kernel of the Markov chain A*^"^ 

We give the proof of Theorem l^[c) when Q!i,...,aAr are distinct. Note that the kernel de- 
pends on Oj in an analytical way, when some Oj are identical, we obtain the kernel by analytic 
continuation. 

Analogous to ( 52 ) , we write the joint probability density formula ( 12 ) of the positive random 
variables A,-"^ into 



p(a«,...,a(^)) 



f n det(0„(AS"\A5."+^)))i:;i,') det(^j^„,(Af )))f,=i, 

n„=ian! Vn=l ,^ J 



where the function (/)„(A-"\ A^"^^'*) {i ^ n + 1) has explicit formula 



x<y, 



and 



{ojv,«jv-i,--.,ajv_j+ 



27ri 



r(2 + 1) ni=i(^ - ctN-k+i 



X 



2 + 1 



-dz. 



(91) 



(92) 



(93) 



Here P(aN,aN-i,-,ciN^_j+i)i^) 3'^^ degree (monic) multiple Laguerre polynomial of the first 

kind, type II (see ( |239[ ) in Appendix [A|). Note that in this section, the functions (!)„,, and 
later are all defined on positive variables. 

We apply Lemma [2] in Section 2.2 with / = [0, 00). Analogous to (63), we have by a similar 
argument that if rii < n2 

.-1 



<n2 •> 



(94) 



where is a large enough contour enclosing all poles an^+i, . . . ,an2, and similar to (64) 



27ri 



27ri 



rZ+l 



X "e -P{a„,...,o„_j_,_i)(a^), 



dz 



Tiz + l) 



-dz. 



j = 0,l,...,n-l, (95) 
j = -l,-2,...,n-A^, 
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where the contour S is the deformed Hankel contour from — oo to — oo that encloses ah possible 
poles of the integrand, z = —1, —2, . . . and z = an+i, . . . , oat counterclockwise. The proof of 
(95) will be given in the end of this subsection. 
Also similar to ( 65 ) , we have 



n 



— aj+i 



Oik 



Note that the vector space Vn is spanned by , x"^ , 

onality (leoj), we have similar to (p| (using ([2371), PMI and p40l)) 



j = 0, l,...,n-2. (96) 
. , Hence by the orthog- 



cI>,"(x)=x-""e-Q(,„,,„_,,„. 



.OLn-j) 



(x) 



X 



27ri /r.r(z + l)nLo(^ 



On-fc 



-dz, (97) 



where Q(a„,a„_i,...,Q„_j.)(a;) is the multiple Laguerre polynomial of the first kind, type I (see 
Appendix [a]) , and the contour Ta encloses the poles a„, . . . ,an-j- 

It is clear again that the Assumption A is satisfied, and the correlation kernel is given by 
(62]) with (/.("i'"2)(x,y) expressed in ([94]), ^l\_k{x) in ([95]) and ^ll_k{y) in 

Next we express the kernel in the double contour integral form. From the formulas (95), 



(97), we compute, analogous to (70) 



n2 
k=l 



dz 



(2«)= h Jr. ^=+'r(«,+i)nf.„,„(---a,)f:^ n£*(«"-°i) 

where Fq, encloses the poles ai, . . . , a„2; and S, the deformed Hankel contour, encloses all the 
poles —1, —2, . . . and F^. Using the telescope trick (71) and (72), we have 

112 



k=l 



{2m) 



X 



-Vr(z + i) /n"ii(^-«i 



niL„2+i(^-"i)^ 



n 



N 

l=ni+l 



{z - ai) 



Note that 



dw 



x-'-^y^T{z + 1) njLn2+i(^ - «j) 



n 



{z-w)r{w + l) Ui=n,+ii- - o^i) 

N 



i=n2+l 



(z - aj) 



n£ni+i(^-aO 

nf=n,+i(^-«j) ^_,,i 
n£ni+i(^-«o 



-x 



-2-1 



F(z + 1) 



r{w + 1) (z - u;) 
F(z + 1) • = 0. 



Thus ( 99 ) can be simplified as 



"2 



fc=i 



(27ri)2 



dw- 



c-^'^y^r{z + i)UT=i(^-'^j) 

{z-w)T{w + i) nrii(^-«o 



(99) 



(100) 



(101) 
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Hence the formulas ( [ToI| ) for Y2=i ^l\-kix)Kl-k^y) and ^ for y) and ^ yield 

the double contour representation of the correlation kernel 



K{ni,x;n2,y) 



2vri Jr. nrin,+i(^-«/) 



(iu'la;>j,ln2>ni 



+ 



-ijy"'r(z+l)njii(^-«i) 
(27ri)2 Jr. "^"^ {z - w)T{w + 1) U?li(.w " «0 



(B dw- 



(102) 



and yield Theorem ^[c) after conjugating out x'^"i /y'^"^, upon proving (95). 



Proof of (jgSj). From (|93j), (l94j) and (|58j), we have 

/•oo 



1 



(27ri)2 



iN 



-dz 



r(u; + 1) nfcLn-j+i(^ - 



(103) 



y 



w+l 



dwdy. 



Here the contour is defined as in ( |94[ ), and the contour S' is a deformed Hankel contour, 
such that they satisfy that for all z G Ta and w G T,' , < ^z. This property will be used 
in (105) to make Jq y^~''"~^dy well defined. Then using the decomposition = — , we 
express ^'"(x) = x""(Part I — Part II), such that 



Part I 



1 



X ^ ^y^ 



(2vrz)2/r.n;i„+i(^-«0 J^' V 

z—l roo 



-dz 



T{W + 1) Uk=n-j+li'W - "fc) 



w+l 



dwdy 



1 



(104) 



where the second identity is the consequence of (pSl), and 



Part II 



1 



X ^ ^y^ 



-dz 



T{w + 1) Uk=n-j+iiw - "fc; 



(2^^)Vr.nL+i(^-«0 
1 



dwdy 



1 

(2^ 



dz 



dwx~''~^T{w + 1) 



n 



N 

k=n—j+ 



l{w - Ok) 



y'-^-'dy (105) 



^T{w + 1) Y\k=n-j-hiiw - ak) 



dz (h dw- 



z — w 



nll„+i(^-ao 



Note that because Pi 



^){y) is a monic polynomial of degree N — n-\-j, the integral 
with respect to y in (104) is F(z+1) times a monic polynomial in z of degree N—n-\-j; because of 
the orthogonal property of -P{ajv,a]v-iv,«n-j+i)(y)' integral vanishes as z = on, • • • , «n-j+i- 
Hence we find that 



Part I 



27ri 



1 / F(z + l)nf=n-i+i(^-afc 



„z+l 



N 



=n+l 



{z - ai) 



dz. 



(106) 
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In Part II, if we integrate z first, by calculation of residues we have 

N / N 



Part n=^(f x-'^-'Tiw + 1) J] _ «,) ^ 



res 



-1 / T{w + l)I[k=n-j+l^'^-(^k 



2m j^, nL+i(^-«o 



(107) 



where we use the identity (res2=^ + X]«ln+i ^^^z=ai){z — w) ^ YliLn+ii^ ~ '^') ""^ = 0- Since as 
integration contours, Fq, + S' = S where S is the contour in (|95|), we have 



Part I - Part II = ^ / r(.; + 1) - 



and we prove (95). □ 



4 Joint distribution of eigenvalues of minor quotients in multi- 
ple JUE 

4.1 Transition probability of the Markov chain A*^"^ 

Analogous to Lemma [T] in Section [2] and Lemma [3] in Section [3| we need the following technical 
lemma. Like in Lemma [3], li,l2, ■ ■ ■ are nondecreasing integers with /„ > n, and an = In — n. 

Lemma 4. Let X' be an In-i x (n — 1) fixed rectangular matrix and Y' he an M' x (n — 1) fixed 
rectangular matrix with M' > n such that R' = {{X')*X' + {Y')*Y'y^{X')*X' has eigenvalues 
jl = {fli, . . . , jln-i) in increasing order. Let X and Y he In ^ n and M' x n rectangular random 
matrices such that 

• The upper-left In-i x (n — 1) hlock of X and the left M' x (n — 1) block ofY are fixed and 
are equal to X' and Y' . 

• Denote the In dimensional column vectorx. := . . . the last column of X, and 
the M' dimensional column vector y := {yi,n, ■ ■ ■ ,yM',n)'^ , the last column of Y . Then 
components of x and y are independent complex random variables in standard normal 
distribution, i.e., p{^Xi^n) = p{^yi,n) = 7\A(0, 1/2) and p{Qxi^n) =p(9yi,n) =7V(0, 1/2). 

• All entries in the lower-left (In — In-i) x (n — 1) block of X are 0. 

Denote the eigenvalues of R = {X*X + Y*Y)-^X*X by X = (Ai, . . . , A„) where < Ai < • • • < 
An < 1- Then the joint distribution of X is 



„-(A)= (^^ + «")' A.(A)niLxAr-(l-A.)^-^'-" 

P^^^^ an\{M' - n)\ A„_i(A) Will f^f-^^l - ^.^/'-n+i ^AdA" ^^^^^ 
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Proof. Instead of R', R, /i and A, we consider T = {{Y')*Y')-'^{X'Y X' that has distinct 
eigenvalues /x = (^i, . . . where < /ii < • • • < fin-i, and T = (Y*Y)~^ X* X that has 

eigenvalues A = (Ai, . . . , A„) where < Ai < • • • < A„,. Below we prove 



^ {M' + a^)\ A„(A)nr=iAr(l + A.)-(^^'+°-+i) 



(110) 



and (109) is a direct consequence of substituting Aj = A.j/(1 — Aj) and = fii / {1 — fxi) into the 
probability measure defined by (110). 

From the assumption of the lemma, {X')*X\ {¥')*¥' and {X')*X' + {¥')*¥' are invertible, 
and from the randomness of x and y, Y*Y and X*X + Y*Y are almost surely invertible. 

Below we assume the invertibility of them. 

By QR decomposition of X', we have Ui G U(/„_i) such that 



X' = Ui 



C 

(i„_i-n+l)x(n-l) 



(111) 



where C is an (n — 1) x (n — 1) upper-triangular matrix. Let U2 = Ui ® G U(/„), we 

have 

C xi' 



X = Uo 



(/„-n+l)x(n-l) X2 



(112) 



where xi = (xi, . . . , Xn-i)"*" and X2 = {xn, ■ ■ ■ , are vectors of dimensions n—1 and Z„— n+1 
respectively, such that if we concatenate them into a /„ dimensional vector x = (xi, . . . ,xi^), 
then C72X = x. We choose an U3 £ U(/„ — In-i) such that 



[73-^x2 = (^,0,...,0)^, where ^ 



and denote C/4 = U2{In-i © t^s), then we have 



In 

E 



^fc,n I ^ ) 



(113) 



X = Ua 



X 

{ln—n)xn 



where X 



C xi 

Olx(n-l) ^ 



(114) 



On the other hand, by the QR decomposition of the M' x (n — 1) matrix Y' , we have C/5 € U(M') 
such that 

YJ{A4' -n+l)x(n-l) J 

where B is an (n — 1) x (n — 1) upper-triangular matrix. Then we have 



Y = U. 



B 



yi 



0{A/'-n+l)x(n-l) y2 



(116) 



where analogous to xi and X2 in (112), yi = {yi, . . . ,yn-i) and y2 = {y ?!;•••) Vm' ) are 
vectors of dimensions n — 1 and M' — n + 1 respectively, and if we concatenate them into 
y = (yi, . . . , yM'), then y = U^^y. We choose an Uq G U(M' — n + 1) such that 



Uq y2 = (ry,0, . . . ,0), where 



M' 



\ k=n 



(117) 
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Analogous to U4, in (114), let = U^{In-i ® U^), we have 

where Y = 



Y = U7 



Y 



(M'~n)xn 



B yi 

Olx(n-l) f] 



(118) 



From ( fTTTj ), ( fini ), ( [TT5| ), ((TTs]), we have 

= C*C, (y')*^' = B*B, T' = {B*B)-^C*C, (119) 

X*X = X*X, Y*Y = Y*Y, T = {Y*Y)-^X*X. (120) 

From the invertibility of {Y')*Y' and Y*Y, we find that B and y are invertible, and then T' 
and T are similar to 



f = {CB-^)*{CB-^), f = {XY-^)*{XY-^] 



(121) 



respectively, which implies the relation between characteristic polynomials, pt>{z) = Pf,{z) and 
Pt{z) =Pf{z). 
We have 

-i_(CB-^ 77-1(^1 -Ci3-iyi) 



XY 



Oi 



in 



-1 



(122) 



•'lx(n-l) 

Since T' and hence T' = {C B^^)* {C B^"^) has eigenvalues /xi, we have the singular 
value decomposition that for Us,Vi G U(n — 1), 

(123) 



CB ^ = C/s-DV]*, where D = diag(^//il, . . . , ^JUn-i)- 
Let C/g = f/s ® /i G U(n) and V2 = Vi ® Ji e U(n), we have 



D ?7 



_V ) where w = [/g-^xi - DF^yi. 



(124) 



We consider the relationship between the characteristic polynomials of T and T' (remem- 



bering pt{z) = Pf{z)), in the same way as (77)-(81), and find 

n-l 



iz- Hi 



-0, 



(125) 



Pt{z) ^ 

PT'{z) 

Pli^ = , _ (^,(A) - a^{^,)) - -(ai(A)ai(/.) + a2(M) - a2(A) - a?(;x)) + 0(1), (126) 



where <ti,o"2 in (|T26|) are defined the same as in (|35|) and (|80|), and Q in (125) are defined 

(127) 



similarly to in (79) 



l-Wjp/ry^ i = l,...,n-l, 



I = n. 



From (125) we find that the eigenvalues of T and T' interlace, that is, A ^ /x. Like (82)-(84) 
in Section |3.1 that there is a homeomorphism between A that satisfies the strict interlacing 
condition ( 82 ) and Ci , • • • , Cn G > and 







res 



in— 1 1 



n-l 



i = 1, . . . , n — 1, 



Cn = o-i(A) - cri(/i) - V Ci = - res 



nLi(^-Afc) _ riiLi^i 



1=1 



(128) 
(129) 
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Hence like (87), we have 



d{Xi, . . . , A„) 



(-1) 



n-l 



A„(A) 



(130) 



We note that the random variables Ci)---)Cn are not independent. However, from the 
definition, components of x and y, and t] and ^ are independent, such that and jji^n are 
in standard complex normal distribution. Thus Wi are independent normal distribution such 
that 3?u)i and Qifj are independent random variables in AA(0, |(1 + distribution, \/2ri is in 
X2{M'-n+i) distribution and \/2^ is in X2(a„+i) distribution. Hence |tfip, . . . , |wn-iP, are 
independent, and their distribution functions are 



1 



1 + /ifc 



e , for i = l,...,n — l, 



1 



(M' - n)! 

By the relation (127) between Q and lifip, . . . , \wn-i 



1 



OLr 



we compute the Jacobian 



9(Ci 



we have that 

p(Ci, 



•,Cn) 



p(Ci, • • • ,Cn,^^ = r-)dr 

fn-\ 



iv 



2\n 



(131) 
(132) 

(133) 



(134) 



M'+a„+l 



poo / ^ \ 

= r'Tn^ir) {'[lp\w,\2{rCi)jPi^2{rCn)dr 

{M' + an)\ C" 

aJ(M' - n)! nr="i'(l + /^.) (i + y^-i ^ , a 

{M'+a^y. nr=iAr"(i+Ai)-(^^'+""+^) 

a„!(M'-n)! nr=iVr(l + /^i)-^''''^"") ' 
Here we use in the last step identity ( |129 ) and the following identity which is a consequence of 



(|128D and ([129D: 

n— 1 ^ 



n-l 



1=1 1=1 



res 



res 



nfc=i(^- Afe) 



res 



n-l, 



=°^(^+i)re=!(^-/^fe 



— iz(z+i)ro-j(z-^fe) 
nr=i(i+A«) 



(135) 



□ 



nr=7(i+/^.) 

Analogous to ( |85| ), we prove plol ) from ( [1301 ) and ([l34]). 

Proof of Theorem^ ( J acobi- Pineiro case). Like in the proofs in the GUE with external source 
and the multiple Laguerre cases, the Markov property of the process relies on that the dis- 
tribution of the eigenvalues of R in Lemma |4] depends on the eigenvalues of i?', but not its 
eigenvectors. By the result of Lemma |4| we have like (46) and (90) that 

p(A(") I A« = i = 



.,n-l)=p^,„-i)(A(")), 



(136) 



where the distribution function p^(n-i) in (136) is defined in (109). Thus we prove the Markov 
property, while ( 109 ) yields ( 10 ) in the Jacobi-Piheiro case, and hence we finish the proof. □ 
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4.2 Correlation kernel of the Markov chain A*^"^ 

Like in Section 3.2 , we prove Theorem [^[d) for distinct ai, . . . , oat, and obtain the general case 
by analytic continuation. 



Analogous to ( 52 ) and ( |9l| , we write the joint probability density formula ( 12 ) of the random 
variables A^"'' whose range is [0, 1] as follows: 

p(a«,...,a(^)) = 

n ^""'^""^^Xr"-""^^' f n det(^4AS'^\A(.-^)))-\) det(^-_^(,(A^)))A..^^^ (^3,) 



n=l 



,n=l 



where the function 0„(A-"\ A^"'"''^'') {i ^ n + 1) has explicit formula 



and 



^JV(^) = nLi (M' -N + j + l + aM-k+i] 



{aN,<y-N-l,--;0iN-3 + l-M'-N) 



27rz 



-z-l 



T{z + 1) ni=i(^ - Qjv-fc+i) 
r(z + j + 2 + M' - A^) 



dz. 



(138) 



(x) 



(139) 



Here Pi 



(ajv,Oiv_i,...,«jv-j+i;^-^'- 



_^)(x) is the j-th degree (monic) Jacobi-Piheiro polynomial of type 
II, and we use the contour integral representation (244) in Appendix [A| and the contour S 
encloses the poles z = —1, —2, . . . , — (j + 1 + M' — N). Note that in this section, the functions 
and later are all defined on variables in (0, 1). 
We apply Lemma [2] in Section 



2.2 



</.("i'"^)(x,y) 



„-i 



with / = [0, 1]. Since our (?i)„(x, y) defined in (138) is very 

1 different, we 1 



similar to the 4>n{x, y) defined in (138), with only the domain different, we have similar to (94) 
that 



~\n2 



(140) 



where the contour Fq, is the same as the in ( 94 ) . Then similar to ( 64 ) and ( 95 ) , 

---ir(z + i)nLi(^-«n-fc+i) 



T{z + j + 2 + M' - n) 
ni=iO' + 1 + M' - n + an-k+i] 



-dz 



2TTi 



{M' -n + j)l 

-n r> 

j-^-ir(z + i) 



(141) 



X X°"fl - x)*^'-"P(,„„..,,„_^^^,;M'-n)(x), 



E r(z + i + 2 + M' - n) nfcIi(-2 - "n+fe 



j = 0,1,. . . ,n - 1, 
■dz, j = -l,-2,...,n- N, 



where the contours S encloses all poles of the integrand. The proof of (141) will be given in 
the end of this subsection. 
Like (96), we have 



n 

.k=j+2 



1 



j = 0, l,...,n-2. 



(142) 
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Then the vector space Vn is also spanned by x"^ , x 



OL2-OLTt 



Hence Hke (66) and 



(97j), by the orthogonahty (|60j) and ( |242| ), ( |243[ ) and ( |245[ ), we have 
(M' -n + j)! 



nLi(i + 1 + M' - n + a„_fc+i) 



(1 — x)" <5(a„,a„_i,...,a„_j;M'-n)(3:^) 



(j + 1 + M' - n + a„_ 



27ri 



x^r(z + i + 1 + M' - n) 

r(z + i)nLo(^-"r^-fe) 



(143) 



where Q(, 



On,an-lv:"n-j 



) (x) is the Jacobi-Piiieiro polynomial of type I (see Appendix IAI) , 



and the contour Fq, is the same as the Tq, in (97). 

It is clear that the Assumption A is satisfied, and the correlation kernel is given by (62) 
with expressed in ([l40j), ^"^.^(x) in ( [l4T| ) and ^'^^.^(y) in pIs] ). To express the 

kernel in the double contour integral form, we write like (70) and (98) that 



k=l 



y-r(z + l) Wtn,+i{w-^i) 



x-+ir(u; + 1) J] 



N 

j=ni+l 



{z - aj) 



"2 



X 



^{M' -k + l + au) 



k=l 

(2^i)2 

n2 



llj=k+ii^ - «j) r{w + M'-k + 1) 

Ulki^-»i) r{z + M'-k + 2) 



dz 



y-'T{z + l)T{w + M' + 1) nL.+i(^-«0 



dw 



2 7r„ x-+ir(u; + l)r(z + M' + 1) J] 



AT 

,i=ni+l 



(Z - Qj) 



fc=i 



(144) 



where the contour S enclose the poles —1, —2, . . . , — M' and the contour Fq (Compare it with 
the E in ([99])). By the telescope trick like ([tT]) and ([72]), writing (z - w){M' - k + 1 + ak) = 
{z - ak){w + M' - k + 1) - {w - ak){z + M' - k + 1), we have 



fc=l 



N 

n 

i=i 



z — a,- 



w — a,- 



AT 

n 



z — a. 



n 



z + M' -l + l 



w-aj -^-^w + M' -l + l 

3=n2+l ^ 1=1 



(145) 



and then 



"2 



k=l 



X "1 y "2 
(27ri)2 



-2-1 «) A^' 

X y TT- w 



na' + fcniil(^-"j) 



dz 



-2-1 w M'-n2 
X y J r 



+ n -ii 



ai 



z — w 



(146) 



Like (|100|), by integrating over Tq, we have 

dz 



-2-1 w M'-n2 

X y -n- w 



dw 

r z — w 



n 



+ /cnj'ii(2-ai 



• dz = 0. 



(147) 
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Hence the formulas ([T46j) for T.T=iKl-ki^)Kl-kiy) and ([TIt]) and ([TIo]) for '"2) (x, y) 
yield the double contour representation of the correlation kernel 



K{ni,x;n2,y) 



-1 

27H 



nrin,+i(^-«/) 



dwlx<yt 



ni<n2 



X "'"2/"' TT- If ■ 



n 



dz <p^^ dw _ ^ ^ ^ U?=i(.w - ai) 



which yields Theorem |i^[d)| after conjugating out /y°"2. 



(148) 



Proof of (141). This proof is parallel to that of (95). From (139), (140) and (58), we have 



-dz 



(149) 



dwdy. 



Here the contour Tq, is defined as in ( 140[ ), and the contour S' encloses the poles — 1, . . . , — (j + 
1 + M' — n) such that 'Rw < ^z for all w £ 12' and z G Fq. Then using the decomposition 
fx ~ lo ~ /o^' express ^'"(x) = x""(Part I — Part II), such that 



Part I 







Uk=r^'iM' -n + j + l + aM-k+i] 



X ^ ^y^ 



-dz 



-W—ll 



T{w + 1) Uk=n-j+iiw - ak) 



(27ri)2 /r. UlLn+ii^ - «0 J^' r{w + j + 2 + M' - n) 



dwdy 



(M' - n + j)! 



2m 



-z-l 



I[f=n+i{z - ai) Jo 



M'-N 



P(a„,at,.i,...,a„-j+i;M'-N)iy)dydz, 

(150) 



where the second identity is the consequence of (139), and 

i X ' ^ ~ y-w-lT-/... , iM-rA^ 



Part II 



v{w + 1) nfcL„-i+i(«' - 



{2-KiY Tr^ nL+i - «/) -^2' F(u; + J + 2 + M' - n) 



1 



dz (h dw 



X 



1 

~ (2^ 
Note that because P/ 



(27ri)2Jr^ Js, r(u; + j + 2 + M'-n) n£n+i(^-«/) -^o 

x"^"iF(w + l) nf=n-j+i(^^-«fc) 1 
^^,—r{w + j + 2 + M'-n) n^In+i(^-«0 



dwdy 

'-""-^dy 



dz d) dw 



(151) 



(ojv,ajv-i,---,On-j+i;M' 



-N)iy) is a monic polynomial of degree N — n + j , 



y'a-y) 



M'-N ] 



(ajv,ajv-i,.--,On-j+i;W 



-N){y)dy 



p{z)r{z + i) 



T{z + j + 2 + M'- n) 



(152) 



where ^(z) is a polynomial of degree M' — n + j, and the residue of p{z)T(z + 1) /F(z + M' — 
n + j + 2) at —{j + 1 + M' — n) is 1. On the other hand, because of the orthogonal property of 
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similar to (|l06l) 



.i,...,a„-j+i;M'-N)iy)i the integral vanishes as z = ajy, . . . ,an-j+i- Hence we find that 



Part I 



27rifr^T{z + j + 2 + M'-n) nL+i(^-«0 



dz. 



(153) 



In Part II, if we integrate z first, by calculation of residues we have like (107) 

/ N 



Part II 



1 



X 



-ji)— 1 



r{w + i) 



N 



2Tii /s/ T(w + j + 2 + M' -n) 



n (^-«'^) E 



res 



'W—l 



r(u; + l) 



1 



27r« Jy,, T{w + j + 2 + M'-n) 



k=n—j+l 
N 



\l=n+l 



n{w - afc) - res 



k=n—j+l 
N 



(^-^)n/ln+i(^-«i). 



-1 r x-'^-^Tiw + 1) nfc=n-j+i - ak) 
2TTi T{w + j + 2 + M'- n) Y{f=n+i {w - ai) 



dw, 



Since as integration contours, + S' = S where S is the contour in (141), we have 

x-^-^r{w + i) Uk=n-j+iiw - (^k) 



Part I - Part II 



2TTi T^r{w + j + 2 + M' -n) Uf=n+iiw - »i) 



-dw. 



and we prove (141). 



(154) 

(155) 

□ 



5 Asymptotic Pearcey process 



For the multiple Laguerre minor process defined by parameters specified in (26), the corre- 
lation kernel for the distribution of the eigenvalues of Sm, Sn2, the minor of X*X, is given by 
(suppose ni,n2 > n) 

-1 / x-'^-^y'" 



K{n,,x;n,,y) = —f _ i^^j ^^l^<.l»i<». 

If,! , x-^-iy'^r(z+l) /z- Inal \" (z- |n6|)"i-" , , 
+ 7 ^ 4> dz (b dw- ' i ^-^^ , (156) 

where Ti^^^jj is a contour enclosing the poles [naJ and [nb\, and Sout is a deformed Hankel 
contour as the S in (21) that encloses Fi^^^j i^^^j. 

For the saddle point analysis later in this section, we need to deform the contour F|^„aj 
into the sum of two disjoint contour Fj^^^j and Fj^^^j that enclose [naJ and [nb\ respectively. 
Then we use the contour Smid instead of Sout) where Smid is, a deformed Hankel contour similar 
to Sout enclosing F|^„^j but not Fj^^^j. 

A simple calculation of residue yields 



dz w dw 



X 



-z-iyWY(^z + 1) f z - [na\ \" (-^ " L»^&J)"'"" 



(2vri)2 Je„.,-e^,, Jr,„^,ur,„,, {z - w)r {w + 1) - [na] J {w - Inb])^^- 



^ dw. (157) 
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From (156) and (157), we have that 



K{ni,x;n2,y) = K^^\ni,x;n2,y) + K^'^\ni, x;n2,y), 



(158) 



where 



K^^'>{ni,x;n2,y) 
,x;n2,y) 



-1 



" \_na\ ^ L J / 

1 



dwln^<n2 =0 \i X <y, 



2TTi 

1 



X -^-^y^ 



(159) 



{w - [nh\ y 



{2m) 



dz 



-dt(;l„i<„2 if a; > y, 

--^-iy"'r(z + 1) [ z - [na\ Y {z - [nb\ )"!"" 



dw 



r[naj ur|^„t,j 



{z - w)T{w + 1) \w - [nb\ J {w- [nb\ )"2- 

(160) 



Although generically na and nb are not both integers, below we use no and nb in place of [noj 
and [nb\ as though they are integers. Readers can verify that it does not affect the asymptotic 
result. 

For the limiting Pearcey kernel to appear, we consider the scaling that ni and n2 are 
2n + 0{n^^'^), and x and y are nc+0{n^^^) where c is defined in (29). After the change of scaling 



(p8|) , we write the kernel of the minor eigenvalue process as (assuming rii = n{2 + cin '^s),n2 
n(2 + cin 2t)) 



K{ni, x; n2,y)dy = — (i^(^) {s, u; t, v) + i^^^) ^. 

(1 + C3n"4n)c-=i" ^(*-*) 



(161) 



where by (158) (noting that we change variables z i— )■ nz,w i— )• nw) 



K'-^\s,u-t,v) 



1 
27ri 



1 — c^n 4 V 



'{b-w)Y^''^^'-^Uwls<t, 



Fa \l — c^n iu J 
if C2 '^""^^(1 + c^n'lu) < C2 '^""^*(1 + csn-it;), 



(1 — Csn 4 7; 
1 - C3n~3n 



(162) 



[c^{b-w)]^'''^^''^Uwl 



s<t, 



ifcs'^" ^^(l+C3n-t^/) >C2'^" ^*(l + C3n-ti 



(27ri)2 



l-can 47;J [c|(5-z)]^i"^^ e^"(^) 1 



dz d) dw- \nz 1 T? I \ f ^ ■ 



(163) 



Here Fa and Ff, are contours enclosing a and b respectively, S is a deformed Hankel contour 
that encloses F^ but not F^, and 

Fn{z) = — \og[cn)nz + logF(n2; + 1) + nlog(z — a) + nlog(6 — z). (164) 

Now we apply the saddle point method to the double contour integral K^'^'^ (s, u; t, v) in ( 163 ). 
Note that by Stirling's formula, for any complex number z satisfying |z| > e and |arg z\ < vr — e. 



Fniz) = nF{z) + ^ log(™z) + O(n-i), 



(165) 



30 



where (taking the principal branch of logarithm) 

F{z) = 2;(log2; — logc — 1) + log(2; — a) + log(2; — b). 
By the definition formulas ( |27[ ) and ( 29 ) , we see that the derivatives 

= log z — logc H ^ h 



F"{z) 
F"'{z) 



z — a 
1 



6' 



z {z- a)2 {z - hf 



+ 



+ 



z^ {z — a)^ {z — b)^ 



(166) 

(167) 
(168) 
(169) 



vanish simultaneously at z = xq. Thus around xq, F{z) = F{xq) + ^F^^'){xo){z — xq)^ + 
0{{z - xo)^). We also have that F^^\xo) < 0. To see this, we denote uq = xo{xo — b) ^ and 
vq = xq{xo — a)~^. Then the fact that xq is the zero of both (168) and (169) implies that 



Let 



ul + v'^ = xo = 2{ul + v^) 
r = uo/vo- 



Then (170) yields 



r(l + r^) 
2(1 + r3)^ 



2(1 + r3)^ 



2\3 



and further 



Xq — a 



Xq 



(1 + r 



2\2 



Vo 2(1 + r3)' 
l + r2)2 



Xo 



1 



_ xo^ 

Uo 

[l + r 



(1 + r 
4(l + r3)2' 

(l + r2)2 



(170) 
(171) 

(172) 



2r(l + r3 



2^2 



(173) 



4(l+r3)2^ 4r(l+r3)2^ 

By the relation a < xo < b, we have uo < and vo > 0, and then r = uo/vq < 0; from the 
relation xo > and its expression in r in (172), we find that r > — 1; and further from the 
property a > and the expression of a in (173), we find that r G (— l,ro ~ —0.657). Now we 
can express F^^\xo) as 



F^^\xo) 



6 



6 



3\4 



16(1 +r3) 



{xo - a)^ {xo - ay 



(1 + r 



(r^ - 3r^ + 8r^ - 3r^ + 1) < 0, (174) 



where the inequality holds for r S ( — 1, —0.4) D (— l,ro) by numerical computation. 

Prom the formula (173), we see that if any one of the three points a,xo,b is fixed, then 

1 + rQ and 



r is determined and so are the other two points. As r 



ro, a 



0, Xo 



'1 — Tq ^)(1 + rp). When r decreases, a,xo,b all increase. If r is close to —1 such that 



-1 + e, then 



Xo 



;e-2 + 0(e"i), xo 



Xo 



;e-' + o(i). 



(175) 



To apply the steepest-descent method to (163), we need the contours Ta^^b, ^ to satisfy the 
following conditions 
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Conditions satisfied by contours Ta,Th,T, 

1. The contours are of the shapes largely depicted in Figure [5} 

2. The point xq is the maximum of ^F{x) for z G S and the minimum for z £Ta and z S Tj,. 
As z moves to oo along any direction of S, ?R.F{x) — )• — oo. 

3. Denoting the contours S, Ff, and Ta near xq by T,^^^, Tf^^ and respectively, such that 
(see Figure [6]) 

X*""" = {z G X I |z-xo| < C3^n"5}, where X stands for S, Fb or Fa, (176) 
we require that the shape of X^^^ is depicted as in Figure |6j 
Below we are going to construct contours as described above. 

Tess ^fiss pess 




XQ+e 4 Co n 3! 



i5n _ -, _ 1 



✓ iTT _ _ 1 



^Ttt _ 1 



Figure 5: The schematic shape of S and 

^a,b = Fa U Ffc. 



Figure 6: The central part of F^*^*^, T,"^^^ and 
F^^^ around xq that are the left, the middle 
and the right respectively. The contours 
extend to the ends having distance n~5 
to Xq. 

The explicit construction of the contours turns out to be tricky. We relegate the construction 
to Appendix [B| and assume the properties above for the contours in the remaining part of this 
section. 



Below we do the steepest-descent analysis of (163) around xq. First we specify the values 
of the constants ci,C2,C3 appearing in psl) as 



1 

3 b — XQ 



C2 = e°-^o, C3 = 6 4,F^'^\xo)i, ci = {b-xofcl. 
For the asymptotic analysis, we write 

K^^\s, U- t, V) = Kjilis, U- 1, V) + i^i'Lnder(^' "5 

where iCess (5, 'u; t, f ) is expressed by the integral formula on the right-hand side of (163) with 



(177) 



(178) 



the contours Fq, F^ and E replaced by F^****, F^**** and T,'^^^ respectively, and K^emainderi^^'^'^^^''^) 
is expressed by the same integral formula with {z, tf) G S x (F^ U F;,) \ T,^^^ x (F^^** U F^*^*^). 
Taking the substitution 



w I— >• xo + Co n iw, z I— >■ xo + Co n ^z, 



(179) 
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such that the contours T,^^^ and T^^, Tf^^ are scaled into E and fa, f ^ respectively (see Figure 
El), we write 



3 \ nxo 

1 — csn 4f \ 



1 — csn 4u 



{b - xo)r^^^ ^'-'^k^l {s, u; t,v)), (180) 



where 



Kislis,u;t,v) 



(27ri) 



dz (p dw 



(1 — Csn 4i;)C3 



(1 — Csn 4u)'^3 



1 -1 -i 

——Co n 4z 



1 

cin 2 ; 



C3 n 4«; J ^ 



1 „-l 



6— xq 3 



cmH gF„(a;o+C3 4«,) ti; — Z 



(181) 



For \w\ < n2o, we have 



1 — CsTl 4,y 



(182) 



and 



—1 — 

c, n 4 10 



1 



6 - Xo 



Co ^n 4tt; 



(183) 



Using approximations (182) and (183) for w € T and z £ T,, and noting that Fn{xo + 
C3 ""^n^^ty) is approximated by nF{xo + Cg "'^n~4ty) as in (165) and 



F{xo + C3 ^n-su-) = F{xo) + ^F^^\xo)ic^\-h)^ + ©((cg in^Jz)^) 

= F(:Eo) + ^z^ + 0(n-Jz5), 
4n 



(184) 



we have 



K^{s,u;t,v) 



1 /• /■ e"^ e"^ 1 g'^Cn^^nto) gCi(n-423) gCi(n-4- ) 

(2^ Jf "^""^ e-^f^ 

...2 „„2 



^ i^z-W (^(n-sn^) ci(„-ii„3) 



(185) 



1 



(27ri) 



dz d) dw- 



e 4 ' 2 



1 



„ I _ 2 I sz 



'*+^-nt w - z 



(l + 0(n-s)). 



Similar to (|180|), we take the change of variable (|179|) and write 

(2) 



remainder 



(s, u; t, v) = c^n-i 



_3 \ nxo 

1 — c^n iv \ 



1 — c^n iu 



:^°(6 - Xo)]'=^"^(^-*)i^i'Lr.der(^, ^)), (186) 



where -f^remamder('^' ^' ^' ^)) defined by a double contour integral formula like Keti{s,u;t,v)) 
in (181), with the same integrand, but the integral is taken on the contour {z,w) G C3n4(S — 



(2), 
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^o) ^ c^ni^Ta^b — xq) \ S X f . Recall that ?R-F{z) attains its unique minimum at xq on S and 
attains its unique maximum at It is not difficult to check that on for (z, w) on the contour 



for K. 



(2) 



) I^Fn{xo+c^ n dominates the other terms 
(2) 



^remainder(^'^;*'^)' ^^e factor e 

of the integrand and make double contour integral -K'remainder(^' ^) vanishing as n — )• oo. 

At last we consider K^^\s, u; t, v) in ( |162 ), the other part of the kernel K(s, u; t, v) in ( 
Note that if s > i, then K^'^\s, u; t, v) vanishes, and if s < t, since ci > and C2 > 1 in d 



161 



177), 



we have C2 "^^^ 



[1 + c^n iu) > 02*^^" ^*(1 + CsTi 41)) for all u,v £ 



and n large enough. 



Below we assume the condition s < t and n is large enough so that 



K^^\s,u;t,v) 




(187) 



Here the contour Ff, like F;, in (162), encloses the pole b. We use a different notation, because 



we want to deform the contour into a square whose left side is through xq, as in Figure [7j so 
that it is of a different shape from the F;, for the asymptotic analysis of K^'^\s,u;t,v). Note 
that Eess is part of F;,. It straightforward to check by explicit calculation that the absolute value 
of the integrand attains its unique maximum on the contour F^ at xq. Then we write similar 



to (178), (180) and (186) 




Figure 7: The square contour FJ around b whose left side is through xq. 



3^0(6 -xo)]'^"^('"*^ 



remainder \ 
]^ _ 1 / 1 — Csn iv * 



X \kgl{s,u;t,v) + K^ 



(188) 



(1) 

remainder ^ 



^{s,u;t,v) 

where (taking the change of variable w ^ xq + c^^n^^w as in ( 179| )) 



s, u; t, v) 



and -^remainder ('^' ^' ^) ^ contour integral with the same integrand as KeU{s, u; t, v) with the 
integral contour changed into Tf, \ T,'^^^. 



2m 



(1 — c^n iv 
1 — c^n~iu 



b-xo 



-1 -i 



(1)/ 



cin ^ 



(189) 
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Using the approximations (182) and (183), we have (if s < t and n is large enough) 



1 

27ri jr 
1 

2TTi 



lOO 2 



(190) 



1 (v^u)^ 1 

e ac-'') (1 + ©(n^J)). 



V'2vr(t - s) 



The estimate of K^^^^^^^^ls, u; t, v)) is similar to -K'^^emamder(^' '^))' ^''^d we can show that it 
vanishes exponentially as n — )■ oo. 

Collecting the results in pHl), (fTTSl), (fTSOl), (fTSSl), plsl), (l90\j and the vanishing properties 



of ■^rei!iainder('5>'";*)''^)> -^'^remainderl'^' ^> and K^si{s , u; t , v) when s > t, we prove that as 
s, u, t, V are fixed and n — >• oo, 

lim (1 + C3ra-3tt)c-^i""^('^~*) {"^/' ^} K{ni,x; n2, j/)dy = {s,u;t,v)dv. (191) 

n^oo Jn{t,V) 

with the conjugation function 



TO/ 



/„(s,n) = (1 + C3n-4n)""o[c-o(^Q _ 5)^] 



-cin2 J 



lim (1 + c^n iu)c — j 77, — )• 00 and n is fixed, 

n— >oo 



Since we have 



and by using the change of variables ( |28| ) reversely 

A(s,7x)=rE"-o((:Eo-6)n)2"-"S 

we prove Theorem |4j 



(192) 



(193) 



(194) 



6 Interpretation in percolation model and Schur process 

In this section we prove Theorem [2] for the three minor processes that we study. We also give 
the proof for the Wishart minor process, for completeness as well as a preparation of the GUE 
with external source case. 

Recall the RSK correspondence that gives a bijection between m x n nonnegative integer 
matrices X = [xij] and a pair of semi-standard tableaux. The shape of the pair of Young 
tableaux, which is a Young diagram, has a bijective relation to the maximal length of / non- 
intersecting up-right paths L^^\m,n) defined in ([6]) with the weight on each site equal to 
the entry Xij of X. A Young diagram can be represented by a sequence of integers in descending 
order, and only the nonzero terms are meaningful. To be concrete, suppose the shape of the 
Young diagram is k^™'") = (k^™'"\ k^2^'^\ . . . ), then 

I 

^kS'"'") =L«(m,n). (195) 
1=1 

Note that since the Young diagram k^™-'") comes from the m x n matrix X by RSK correspon- 
dence, K^™"'"^ = for all i > min(m,n). 
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Below we fix m and assume 1 < n < m. By definition, k^'"'") has infinitely many component, 
but since the components kI!"'"^ = for k > n, we abuse the notation and consider as an 



■variate vector, as a discrete counterpart of the /i^™-'") defined in ([s]). From the relation (195) 



to the last-passage paths, we have the interlacing property k^"^'" ^) ^ k^™'""-*, that is, (note that 
as i increases, our k^™"'"^ are decreasing but /i^™'"'* are increasing) 

U ^ K„ ^ '^n-l — — — — — ^2 — '^l — '^1 ^ OO. 

(196) 

In other words, {k^™'"^ | 1 < A; < n < A^} form a Gelfand-Tsetlin pattern. 

Let si, . . . , Sm', h, . . . ,tn be positive parameters such that Sitj < 1 for alH = 1, . . . , m, j = 
1, . . . ,n and let w{i,j) be an independent random variable in geometric distribution with pa- 
rameter Sitj, so that 

F{w{i,j) = k) = {l- s^tj){sitjf. (197) 



If the entries Xij of the mxn matrix X are random variables w{i,j) defined in (197), the shape 
of the corresponding pair of the Young tableaux, which is a random Young diagram, follows the 
Schur measure [22], i.e., the probability that the shape is k is 

m n 

(ti,...,t„). (198) 

i=ij=i 

Equivalently, if the distribution of the weight w{i,j) at each site of the mxn square 

lattice is defined by (197), the maximal length of I non-intersecting up-right paths L^^\m,n) 



becomes a random variable and is determined by (198) via (195). Note that the maximal length 
L(')(m,n) is well defined if the weight Xij are real numbers instead of integers. Then we can 
define nf^'"^ conversely by ( |195 ). In that case = {k,^^'^\ k^^'"\ . . .) is a sequence of 



real numbers in descending order, but may not be interpreted as a Young diagram (cf. 
defined in ([s]) and Remark [7|) 

If we fix m and let n vary between 1 and A^, we have the joint probability distribution of 
. . . , as follows. 

Proposition 1 (pOj). In the m x N rectangular lattice \ i = 1, . . . ,m, j = 1,. . . ,N} 



where N < m, let the weight at each site be independent random variables given by (197). If 
n runs from 1 to N, and is regarded as the discrete time, the random Young diagram k^™-'") 
evolves as an inhomogeneous Markov chain with transition probability 

P„_l,„(At('"'"-^),K('"'")) = 

S,,^,„,{si,...,Sm) ,EL,>c^'"'-Eri^«^-"-^\ , _ ^ .-iqoA 

~Jl 0^{m,n-l)(Sl, . . . , Smj 

and the initial state = (0). 

In the special case that qj = 1 for all j, Proposition [T] is proved in [2T]. See also [12]. 
Choosing special values oi pi,qj and taking limit. Proposition [T] yields the four cases of 
Theorem [2] Below we show the detail of taking limits case by case. 

Wishart Let the m = M in Proposition [T] Set 

Pi = 1, Qj = 1 + Oj/L, for i = 1, . . . ,M, j = 1, . . . ,N, where aj < 0, (200) 
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we have 



1- w{i,j) 



(201) 



L— >oo L 

where u{i,j) are independent random variables in exponential distribution with parameter —aj, 
as in the continuous percolation model in the Wishart case of Theorem [2] 

Note that in this special case ( 199 ) can be simplified by the evaluation of Schur polynomials 
[161 Section 6.1, Exercise 6] 




Am(1 - Ki,2 - K,2, . . ■ ,M - km) 
Am(1,2,...,M) 



(202) 



if KM 



[Lx] 



+1 

{M,n) 



B If 

(M,n) 



I If we take the limit L — t- oo in (199) and consider k. 



{M,n) 



KM+2 = • • • 

for fixed values of x\'"'"' , we find that the transition probability (199) in Proposition 



[l| satisfies (LL:e(A'^'")J = ([Lxi*''")] , . . . , [Lxi^''''^ \)) 



lim L"P„_i,„([Lx(*^'"-i)j, L^x(*'^'")j) 



(-an) 



M 



Anix. 



{M,n) 



{M,n)s 



(M,n) 



(M-n)!A„_,(4iY''^ 



, x^ 



(M,n-1) 



and 



lim LP{lLx[''''^ \) = (-ai)^'^((M- l)!)-i(xi''''V^-ie''i^ 



L— >oo 



(M,l) 

^ 1„(M,1) 



>0' 



Now if we consider random variables ^^'"^ defined by ^ as 



(M,n) 



(M,n) 
LA ' 



for 1 < i < n < iV, 



(203) 
(204) 

(205) 



then as L — )• oo, converges in distribution to /ij-*^'"^ that is defined in the continuous per- 



colation model by ([8]) where u{i,j) are independent random variables in exponential distribution 
with parameter —aj. Then constitute an inhomogeneous Markov chain as n runs from 1 to 

N, and the transition probability P„_i,„(//(^''^-i), is lim^^oo L'^Pn-i.nd^Ai^^'""^^] , L-^^/x^^-^'^^J) 
given in (203), and the distribution of ^J.{^^'^^ is liniL^oo LP {[^[^^'^^ \) given in ( |204 ). We see 
that the transition probability of coincide with the transitional probability (10) in the 

Wishart case, and the probability distribution of /i^*^'^'' is the same as the pi{fi\^^'^^) in Table 
[TJ Thus we prove the Wishart case of Theorem [2| 

Remark 7. The quantity fJ-\^'"'^ that come from the limit RSK correspondence can also be 
associated to the continuous RSK correspondence. See [19], [20] and [HI Appendix A] for more 
details. 



GUE with external source The classical transition between the Laguerre and Gaussian 
weights gives that if n is finite and x = M + \/M^, 



lim e^-'^^^M-^+'^e-(i-'^/^^^)^x^'^-" = e-V+<. (206) 

Suppose in the M x N rectangular lattice {(i, j) | i = 1, . . . , M, j = 1, . . . , N} [N < M), the 
weight at each site is an independent random variable u{i,j) in exponential distribution with 



Vhere Am{xi,. . . ,xm) = n"i H^ii+ife " ^i) 
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parameter 1 - ajj^fM. Let /x^^'^), . . . be defined by ^ and for 1 < z < n < TV, we 

define the random variables = (zyj^'"-*, . . . , v^^'"^^^ by 



;,(*^'-)=M + VMz.f''"\ (207) 

such that X^Li i^f^'"'' = M~^/'^{L^^\M,n) - IM). From the property of we have that 

as n runs from 1 to A^, the sequences i/^^'^^ = . . . , u^^'^^) evolves as an inhomogeneous 

Markov chain. Moreover, as M — )■ oo, i/^-^'"') have a joint limit distribution. To compute the 
limit, consider the transition probability Pn-i,n{fJ'^^^''^~^\ fJ'^'^'^'^^) (10) in the Wishart case of 



the last-passage percolation model with parameters specified by pi = 0,qj = 1 — aj/\/M, we 
have, substituting ^(*^'") into by ([207| and assuming = for ah M, 



, (n-1) , (n-l),2 l^(n-l)^^(n) (^iUSj 

27r A 1 f 7/'^ lA ') 1 , (n-l) 



nr=i 

and 



lim VMPifir'^) = (2^)-i/2e-'^?/2e-iW -"^)\ (209) 

A/— !>oo 

A simple scaling argument yields that as M — t- oo, the limiting transition probability from 
j^{M,n-i) j,(M,n) ^las the density function limM-^oo M"/2p„_i_„(/i(A'^'"-i), ^(^'-f.'^)), and the 

limiting distribution of is limjv/-5>oo \/MP(/i^*^'^''). We see that the limiting transition 

probability of coincide with the transitional probability (10) in the GUE with external 



source minor process and the limiting distribution of z^|^^'^^ is the same as the pi{z) given in 
Table [T| and so we prove the GUE with external source case of Theorem [2] 

Jacobi-Pineiro Let the m = M' in Proposition [!} Set 

Pi = t'-\ qj=r^+\ forz = l,...,M', j = l,...,iV, (210) 

where t £ (0, 1) and aj > 0. By the Jacobi-Trudi formula of Schur polynomials, for any Young 
diagram k = (ki, K2, ■ ■ ■), if km'+i = km'+2 = • • • = 0, then |16l Section 6.1, Exercise 5] 

_ det(t(-i)K-+^^'-^-))i^L, _ AM'{t^^^'''-\t^^+'''-',...,t^^") ^ , 
' ' ^ det(t(-i)(*^'-i))MLi AM'(^*''-^^^'-^ •••,!) ' ^ ^ 



Thus with our special choice oi pi,qj (210), the transition probability (199) becomes 



P 1 (JM\n-l) (M',n)^ _ TJ(. ^M'[t , ■ ■ ■ ,t ) 

l-n-lMI^ ,K ) - Lly^ ^ ' . , (M',n-l)^,,, (M',n-1)_^.,, „ ^(M,n-1) 

i=l ^M'{t''^ >+M'-1^^k\ ) 

I iiNz-V^n (Af'.n) (M',n— 1)^ 

Xt{"n + 1)(E.=1-1 ^-E.=l«l )l,,M'.„-l)^,(M',„), (212) 



where we assume '^^ = k"^^ = for i > n,j > n—1. Note that k-^^ '"'^ +M' —i = M' — 1 

JM',n-l) 



(resp. K^- ' + M' — j = M' — j) if i > n (resp. j > n — 1). 
Now we assume 

t = e-i/^, (213) 
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Similar to (203), if we take the limit L — )• oo in (199) and consider '"^ = \_Lx'^^ '""^J, using 
the identity 



(Af',n) , , . (M',n) 

lim = e-^« , for 1 < i < n < A^, 



(214) 



we find that the transition probability (199) in Proposition [T] (( [212 ) in our special case) satisfies 



lim L"P„_i,„([Lx(*^''"-i)j, Lix(*^''") 



^ M' 



{m' ,n) (M',n)^ ^„ , (A/',n) ^ , ^ (M',n) 



{M',n-1) (M',n-1) 

An-i(e ^1 ,...,e ^"-1 )ni=i (1-e 



(A/',n-l) ^ , , . , ,-. (M',n-1) 

-xl ' jAf'-n+lg-(On + l)x> ' 

Also we have 

M' 

lim LP(Lxi*'' ^)j) = ((M' - l)!)-i TT(«i + - e^'" ■^y^'-ie-("i+i)4' 



i=l 



On the other hand, since t depends on L by (213) similar to (201), 



(215) 



(216) 



hm — - — = u{i,j), 

L— >oo Li 



(217) 



by (210) is a random variable in exponential distribution with parameter aj + i. Then like in 



the proof of the Wishart case of Theorem |2 , we consider random variables '"^ defined by 



(M',n) 

as 



(M',n) r (A^'=") f 1 ^ ■ ^ ^ 7\r 
= -ty/ij^ j , tor 1 < ? < n < A'. 



(218) 



As -L — 7- oo, '"^ converges in distribution to /i-*^ '"^ that is defined in the continuous perco- 
lation model by ([s]) where u{i,j) are independent random variables in exponential distribution 
with parameter aj+i. Then /i^*^ '"^ constitute an inhomogeneous Markov chain as n runs from 1 
to N, and the transition probability P„_i,„(^(*^>-i), /i^*^''") is lim^^oo i:''P„-i,„(L^^^^^''""^^J , L^/^^*^''''^] ) 
given in (215), and the distribution of '^^ is liniL^ao LP {[iJ.[^^ '^^ \) given in ( |216[ ). To prove 
the Jacobi-Piiieiro case of Theorem [2j we need one more step of change of variables, such that 



(M',n) -a 



(M' ,n) 



for 1 < i < n < iV. 
Noting that constitute an inhomogeneous Markov chain and 



(219) 



n—l,n 



(M',n-1) ,.(M',r^)^ 



i=l 



(M',n) 



p(^(A^M))^p(^(MM))^(AfM)^ 



(220) 
(221) 



we check that the transition probability of z^(a^''") coincides with the transition probability 

(M',1) 



(jlOj) in the Jacobi-Pirieiro case, and the probability distribution of is the same as the 

pi(^^^^ '^'*) in Table [l| thus finishing the proof in that case. 
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Multiple Laguerre limit The classical transition between the Jacobi and Laguerre weights 
gives that if ^ is finite and x = S,/M' , 



lim (M')°x"(l-x 



,M' 



(222) 



Then similar to the derivation of the GUE with external source case from the Wishart case, we 
find the proof of the multiple Laguerre case of Theorem [2] from that of the Jacobi-Piiieiro case. 

Suppose in the M' x N rectangular lattice \ i = 1, . . . , M' , j = 1, . . . , N} {N < M), 

the weight at each site is an independent random variable u{i,j) in exponential distribution 
with parameter aj + i. Let /U^^^ . . . , ^^^^ be defined by ([S]) and for 1 < i < n < A^, we 
define the random variables = (p^*^ '"'^) and zy(*^''") = {v^^ . . . , '"'^) 

by 



-(M',n) 



(M' ,n) 



{M',n) 



m'd} 



(M',n) 



(223) 



such that X]i=i ~ —L^^\M',n) + llogM'. From the property of fj,^'^^ we have 

that as n runs from 1 to N, the sequences = (i/j*''^ . . . , z>^'^^ '"^) and i/^*^''*^) = 

(t'n^ • ■ ■ , i^n*^ '"■*) evolves as an inhomogeneous Markov chain. Moreover, as M' — )• oo, 
have a joint limit distribution. To compute the joint limit distribution of 

jy{M\n)^ we 

note that from the result obtained above in the Jacobi-Piiieiro limit, the distribution function 
of is given by (10) in the Jacobi-Pirieiro case, with ^(^'''^''") substituted by 

^{M',n-i)^^{M',n)_ Similar to (Iffl, if we fix z/(*^''") = j/W for all M' and take M' ->■ oo, then 



lim (M')-"Pn™i,„(p(^^'"-i\z:^(*^'")) 



in). 



OLr. 



(n-l) 



Yyl^(n-l)-<^(n) (224) 



and similar to (209) 



(225) 



Compare the limiting transition probability (224) with the transition probability (10) in the 

multiple Laguerre minor process and limA/'^oo(M')-ip(i>P''''^) with pi(z) given in Table [ll 
we prove the multiple Laguerre case of Theorem [2] 



A Contour Integral Representations of multiple Laguerre poly- 
nomials of the first kind, and Jacobi-Piiieiro Polynomials 

The multiple orthogonal polynomials are multiple weight generalizations of orthogonal polyno- 
mials, such that the single weight that defines the orthogonality is replaced by a sequence of 
weights. Among the orthogonal polynomials, the Hermite, Laguerre and Jacobi polynomials 
are the most well known, and are called very classical polynomials. The three kinds of very 
classical orthogonal polynomials have various multiple weight generalizations. With our scope 
being restricted to the so called AT systems, there are four corresponding very classical multi- 
ple orthogonal polynomials. To be concrete, the Hermite polynomials give rise to the multiple 
Hermite polynomials, the Laguerre polynomials have two generalizations, namely the multiple 
Laguerre polynomials of the first and second kinds, and the multiple weight generalization of 
the Jacobi polynomials are called Jacobi-Piheiro polynomials. The algebraic properties of these 
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four kinds of multiple orthogonal polynomials are summarized in [23] , and the relation among 
them is illustrated in the diagram in the beginning of Section 3 of [23j. 

The multiple Hermite polynomials and the multiple Laguerre polynomials of the first kind 
(simply called the multiple Laguerre polynomials in some random matrix theory literature) are 
closely related to the random matrix models called GUE with external source and LUE with ex- 
ternal source (better known as complex Wishart ensemble) respectively [B] . The random matrix 
models corresponding to the other two kinds of very classical multiple orthogonal polynomials, 
as far as we know, are first studied in this paper. For the sake of asymptotic analysis, it is 
desired to have contour integral formulas for these multiple orthogonal polynomials. To our 
limited knowledge, these contour integral formulas have not been written down explicitly in 
literature. We state and prove the formulas in this appendix. 

For comparison and reference, we first state the contour integral formulas of multiple Her- 
mite polynomials and multiple Laguerre polynomials of the second kind, both obtained in [6]. 
Unlike orthogonal polynomials with respect to a single weight, there are type I and type II 
multiple orthogonal polynomials with respect to a sequence of multiple weights. First we con- 
sider multiple Hermite polynomials. Let ai,...,an be a sequence or distinct real constants. 
The n-th multiple Hermite polynomial of type II, denoted as P{a-i,...,a„)i^)i is an n-th degree 
monic polynomial satisfying the orthogonality condition with respect to the multiple weights 

2 2 

/oo 2 
Pia„...,a„)(.x)e-'^+'"^''dx = 0, k = l,...,n. (226) 
-oo 



The (n — l)-th multiple Hermite polynomial of type I, denoted as Q(ai,....a„){x)^ is not a poly- 
pi ' ^2 ^2 

nomial (see Remark^below) but the linear combination of the weights e"^"*""^!^, . . . , 
satisfying the orthogonality condition 

/oo 
Q(au...,a^){x)x^dx = 0, k = 0,...,n- 2, (227) 
-oo 

Q{a„...,a^){x)P(a„...,a„^,){x)dx = 1. (228) 



' — oo 

fOO 



There are contour integral formulas for the multiple Hermite polynomials of both types I and 
II. 

Proposition 2 ([6j Theorems 2.1 and 2.3). 

where the contour C -|- iM t is the upward vertical line {z = C + iy \ y £ M}, and 



If " 1 

Q(„, „ ){x) = —= <t e-i(*-^)' TT dt, (230) 

where Ta encloses all poles ai, . . . counterclockwise. 

Remark 8. Our definition of the multiple orthogonal polynomials of type I is an abuse of lan- 
guage, but it is essentially equivalent to the standard definition. (Our Q[ax,...,an){x) corresponds 
to the Qfi{x) in Formula (1.6) in [6].) This remark also applies to other multiple orthogonal 
polynomials of type I defined below. 
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Remark 9. Our Proposition [2] is only the generic case of the theorems in [6] in the sense that our 
-P(ai,...,a„) and Q{ai,...,an) ^rc the Pfi and with n = (1, . . . , 1) in [6]. to recover the degenerate 
cases, it suffices to take the hmit that some become identical. 

The multiple Laguerre polynomials are defined similarly. Let p > —1 be a real constant and 
oi, . . . , a„ < be a sequence of distinct negative real constants. The n-th multiple Laguerre 
polynomial of the second kind, type II with respect to the multiple weights 
denoted as P{ai,...,an;p)i^)^ n-th degree monic polynomial satisfying the orthogonality con- 
dition with respect to the multiple weights 

/•CO 

/ Pia„...,a,.;p){x)xPe'''''dx = 0, k = l,...,n. (231) 
Jo 

The (n— l)-th multiple Laguerre polynomial of the second kind, type I, denoted as Q(ai,...,an;p)ix): 
is the linear combination of x^e"^!^, . . . jX^e"""^ satisfying the orthogonality condition 

Q{a,,...,a^;p)ix)x^dx = 0, /c = 0, . . . , n - 2, (232) 

/•oo 

/ Q{ai,...,ar,;p){x)P(ai,...,ar,.i;p){x)dx = l. (233) 

Jo 

There are contour integral formulas for the multiple Laguerre polynomials of the second kind, 
both types I and II. 

Proposition 3 (^6] Theorems 3.1 and 3.2). Suppose p is a nonnegative integer. For x G (0,cxd), 

Pia.,...,a„;p){x) = ^(n + ^ + 1) ^ I ^n-p^l "Q^^ _ (234) 

1 ik=l "-k 

where Ti is a contour enclosing counterclockwise, and 

T-rn P r " 1 

Q(a, a ■v)ix) = (-l)^i#^— (t e^'*r+f-i TT ^—dt, (235) 

where Ta is a contour enclosing poles ai, . . . , a„ counterclockwise. 

Note that our contour integral formulas (234) and ( 235[ ) are different from [6^, Formulas (3.5) 
and (3.10)] by change of variables. An analogue of Remark [o] applies to Proposition [sj 

Remark 10. Although the multiple Laguerre polynomials are well defined for all real p > —1, 
the statement and the proof given in [6] of the proposition are only valid for integer p. When p 
is not an integer, the expression (234) is valid if we let S be the deformed Hankel contour, but 
there is no simple way to generalize the contour integral formula in (235) for noninteger p. 

Now we state the result for multiple Laguerre polynomials of the first kind. Let ai, . . . , > 
— 1 be a sequence of distinct real constants. The n-th multiple Laguerre polynomial of the 
first kind, type II, denoted as P(ai,...,a„)ix), is an n-th degree monic polynomial satisfying 
orthogonality condition 

/•oo 

/ i'K,...,a„)(x)x"'=e-"dx = 0, k = l,...,n. (236) 
Jo 

The (n — l)-th multiple Laguerre polynomial of the first kind of type I is the linear combination 
of x^^e"^', . . . , x""e~^ satisfying the orthogonality conditions 

/•oo 

/ Q{a„...,a,.)ix)x''dx = 0, k = 0,...,n- 2, (237) 
Jo 

Q{au...,au)i^)P{a^,...,a„-i){x)dx = 1. (238) 
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There are contour integral formulas for the multiple Laguerre polynomials of the first kind, 
both types I and II. 

Theorem 5. For x £ (0,oo), 

where the contour S is the deformed Hankel contour, and 



27rz /r^ r(z + 1) nfc=i(^ - Ofc) 

where the contour Ta encloses the poles ai, . . . , a„ counterclockwise. 

At last we state the contour integral formula for Jacobi-Piiieiro polynomials. Let /3 > — 1 
be a real constant, and ai,a2, ■ ■ ■ ,an > —1 be distinct real constants. The n-th Jacobi- 
Piheiro polynomial of type II P(ai,...,an;i3)i^) is an n-th degree monic polynomial satisfying the 
orthogonality condition 

^K,...,a„;/3)(x)x"Hl-a^)^rfa; = 0, k = l,...,n. (241) 

The (n — l)-th Jacobi-Piheiro polynomial of type I, denoted as Q{ai,...,an;i3){x)i is the linear 
combination of x°'=(l — x)^, k = 1, . . . ,n, that satisfies the orthogonal conditions 

/ Q{a^,...,a„M^)x''dx = 0, k = 0,...,n- 2, (242) 
Jo 

1 

Q{a-,,...,a„;l3)ix)P^au...,c.„_ul3)ix)dx = 1. (243) 

There are contour integral formulas for the Jacobi-Piheiro polynomials of both types I and 

II. 

Theorem 6. For x G (0, 1), 

r(n+l + /3) jl-x)'^ /- ^-^-ir(z + l)nLi(^-«fe) ,.. .244^ 

where the contour S is the deformed Hankel contour, and 

"^(--■--^^^"^ - ^ 1 r(z+i)nLi(--«.) ' ^ ^ 

where the contour encloses the poles ai, . . . , a„ counterclockwise, but does not enclose other 
poles, if there are any. 

We note that if f3 £ the contour S can be a closed contour since the integrand has only 
finitely many poles, and Ta can be any contour enclosing ai, . . . since the integrand has no 
other poles. 

In the case n = I and /3 + ai = —1, on the right-hand side of (245) the constant factor 
vanishes while the contour integral blows up. However, (245) still holds via VHopital's rule as 
ai — ;> /3 - 1. 

Before the proof of Theorems [5] and [6j we recall the following results about gamma function. 
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Lemma 5. Let C and e be positive constants, and suppose < 0, < C. 

(a) If \z + k\ > e for all k = 0,1,2, ... , we have 

T{z) =x-^-ie^O(l), 
and the 0(1) factor is bounded uniformly in z. 

(b) If a,b gR. and \z + a + k\ > e, \z + b + k\ > e for all k = 0,1,2, ... , we have 

r{z + b) ^^>' 

and the 0{1) factor is bounded uniformly in z. 
Proof. Applying the Euler reflection formula (see [ii Theorem 1.2.1]) 

T{z) 



(246) 



(247) 



sin(7r2;)r(l — z) ' 



(248) 



we can derive the desired asymptotics of T{z) (resp. T{z + a) and V{z + b)) from those of 
r(l - z) (resp. T{l-z- a) and r(l - z - b)). Note that "Rz <0 implies 5R(1 - z) > 1. Apply 
the asymptotic formula Corollary 1.4.3]) 



V{w) ~ V27rw"'"2e""', for |argu;| < vr - 5, where 5 G (0,7r). (249) 
we obtain the asymptotics of r(l — z), and hence prove ^[a) The proof of ^[b)| is similar. □ 

In the proof of the theorems, we use the contour to realize the deformed Hankel contour 
S, defined as 



^a,b = ^a,b U ^a,b U ^a,b, ^here ( 



' fj = {z = t — ai \ —oc < t < b}, 

= {z = b + it \ —a <t< a}, 
^ Y,^ fj = {z = ai — t \ —a < t < oo}. 



(250) 



Proof of Theorem [s] First we prove ( 240 ) . A simple residue calculation shows that 

res — , , J!^„ — ; -dz (251) 



ak r(z + i)nLi(^-«fc) 



is a constant multiple of Thus the right-hand side of (240) is a linear combination of 

^Qfcg-x^ A; = 1, . . . ,n. 

On the other hand, for any m G Z+, 



-dz x^dx 



Jr, r(z + 1) nLi(^ - 



r(z + m + 1) 





1 

27d 

1 

2vri Jr, r(z + 1) U]=iiz - aj 

-— (k Fm{z)dz, where Fm[x) = ™ — 

27rUr, Uj=iiz-oij 



r(^ + i)nLi(^-«. 



-dz 



-dz 



(252) 
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For m = 0, 1, ... ,n — 2, Fm{z) = 0{z'^~^) = 0{z~'^) as z — )• oo. Deforming Fq into a large 
enough circle entered at the origin, we see that 2^ §^ Fm{z)dz vanishes, and prove that the 
right-hand side of (240) satisfies (237). If m = n — 1, we have Fm{z) = z'^ + 0{z-'^). With 
Fq deformed into a large circle centered at the origin, we see that ^ Fm{z)dz = 1. This 
result, together with the vanishing result obtained above for m = 0, . . . , n — 2, yields that the 
right-hand side of (|240l) satisfies (|238|). 



Next we prove (239). To show that the contour integral over S in (239) is well defined, we 
apply Lemma c|a) , and find that for x ^ Q, the integrand of the contour integral in (239) decays 
rapidly as z moves along the deformed Hankel contour and 'Rz — )• —00. Thus ..,«„) (x) is 
pointwise defined for x 7^ 0. 

Since F(x) has poles z = 0,-1,-2,..., the poles of the integrand in ( |239 ) are z = 



-1,-2,, 



At the pole z = —j, we denote the residue 

r(^+i)nLi(^-«^) _ (-irnLi(j+«fc) 



res 



,.2+1 



(i - 1)! 



(253) 



Thus we can formally apply the residue theorem, and express the contour integral over S in 
(239) as ^ji ^ power series in x. To make the argument rigorous, we write 



1 / r(z + i)nLi(^-«fc) 



r.2+1 



dz 



l. — k—7j j — 1 



we have that as /c — )• 00, 



where is defined in (250). li x £ (0,oo), by Lemma 

the contour integral over Ti^ vanishes for any x 7^ 0. From ( |254 ), we find that X^j^i 

converges for any 2; G C. Thus the contour integral over S in (239) is a power series in x for all 
X G (0, 00). 

To show that the right-hand side of (239 ) is a polynomial in x of degree n — 1, we apply the 
identity that for k £ Z, 



1 

2TTi 



T{z + k + l] 

rpZ + l 



dz = X e 



(255) 



To prove (255), we note that the contour integral in (255) is similar in form to that in (239). 
By arguments like above, we express the contour integral in (255) in power series and have 

r(z + k + i 



1 

27ri ./2 



-dz 



E, 

j=o' 



res 

-j-k-i 



T{z + k + l) 



E 

j=0 



{-iyx^+'' 



x'^e' 



(256) 



From the recurrence formula of the gamma function, we find that there are coefficients 
Co, . . . , Cn such that 



Thus 



r(^ + i)nLi(^-«fe) 

0-2+1 



k=0 



F(z + A: + 1) 



X 



2+1 



where c„ 



1. 



e 

27ri 



r(^ + i)n-=i(^-«;^: 



r.2+1 



dz = CfcX^ 



A:=0 



is a monic polynomial of degree n. 



(257) 



(258) 
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We need to show that the right-hand side of ( 239 ) satisfies the orthogonahty relation ( 236 ) . 
Given any M > 0, we have for any m = 1, . . . ,n 



M 







1 / r(z + i)nLi(^-«fc) 



27ri 



rZ + l 



1 
27ri 



x^-^^^^dx J r{z + 1) JJ(z - ak)dz 
^ k=i 

- I W'-^-^Tiz + 1) TT (z- ak)dz 
"i It 

•'^ k=l,...,n 



(259) 



27ri 



n-l ^, 



— (p M°'"-^r(z + j)dz 

^ 27rz is 



where the constants c'^, . . . , c'„_^ are similar to the constants ci, . . . , c„ in (257) such that —Vi^z-^- 



1)11/0=1,.. .,n(2 - Ofc) = E?=i c'r(z + j). By (p55|), we have 



k^m 



1 

27ri 



M"'"-^r(z + j) =M"'"+^e 



Substituting (260) into (259), we verify the orthogonality condition by 

r-M 



lim 

M-5-00 







r(^ + i)nLi(^-«fe) 



dz x'^'^dx = 0. 



(260) 



(261) 



Proof of Theorem [6| First we prove (245). Without loss of generality, we assume P + ai ^ 1 
if n = 1. A simple residue calculation shows that 



res 

z=ak 



x^r(z + n + /3) 

r(^ + i)nLi(^-«fc) 



(262) 



is a constant multiple of x"''. Thus the right-hand side of (|245|) is the linear combination of 

x'^ 



{l-xf, k = l,...,n. 
To prove the orthogonality condition (242), we note that for any m = 0, 1, . . . , n — 2, 

f{l-xf r x^T{z + n + f3) 



dz x^dx 



2-Ki JT^^{z + l)YYk=i{z- ak) 

r(z + n + /3) 



r(^ + i)nLi(^-«fc 

r(/3 + l) r T{z + m + l)T{z + n + P) dz 
27ri Jr. r(z + m + 2 + /3)r(z + 1) OLi " 



(263) 



r(/3 + i) 

27rz 



Gm{z)dz, where G.m{z) 



nr=i(^+j)nr=>;+2(^+^+/3) 
nLi(^-«fc) 



Similar to (252), for m = 0, 1, . . . , n — 2, Gm{z) = 0{z^ ") = 0{z ) as z — )• oo. By deforming 
into a large circle, we have that 2^ Gm{z)dz vanishes, and prove that the right-hand 
side of (I245I) satisfies (IMl. 
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To verify the orthogonality condition (243), we have similar to (263) 
1 







dz 1 ^dx 



[l-x^ r x'T{z + n + f3) 
2m /r„r(z + l)nLi(^-«fe) 

= r(/3 + l)* Gn-i{z)dz, where G„_i(z) 



ni=i(^ + ^') . (264) 

(z + n + /3)nLi(^-"fc) 



If we replace the contour into a large circle F, by the estimate 



1 
27ri 



= 0, 



0{z ^) we have 
(265) 



Note that the contour F encloses one more pole z = —{n + (3) than the contour Fq. We have 

n-=i(^+i) 



by (265) that 



2m 



Gn-idz 



res 



n-=i(j+/3) 



(n+/3) (^ + n + /3)nLi(^-«fc) riLil'^ + Z^ + afe)' 



(266) 



By (264), (266) and the vanishing of ^ Gm(z)dz for m = 0, 1, . . . , n — 1, we verify that the 
right-hand side of (245) satisfies ( 243[ ). 

Next we prove (244). To show that the contour integral over E in (244) is well defined, we 
apply Lemma [£|[b)| and find that for x 7^ 0, we see that the integrand in ( |244 ) decays rapidly as 
z moves along the deformed Hankel contour and — )• —00. Thus Q[ai,...,a„;i3){^) is pointwise 
defined for x / 0. 

Since T(z) has no zero and has poles at 2: = 0, —1, —2, . . . , the poles of the integrand in 
(244) are z = —1, —2, .... At the pole z = —j, we denote the residue 



X 



-z-l 



res 

z=-j 



r(^ + i)nLi(^-«fc) 



res 

z=-j 



F(z + n + 2 + /3) 

x— iF(z + i + l)nLi(^-«fc) 
T{z + n + 2 + I3)lli=i{z + l) 



(267) 



^ (-1)"-^' m=iU + (^k)ULiin-l + 2 + /3) 
F(n + 2 + /3) (j-1)! 

Thus we can formally apply the residue theorem, and express the contour integral over S in 
(244) as Yl'j^i -^j' ^ power series in x. To make the argument rigorous, we mimick the argument 
for the contour integral in ( |239 ) and write 



1 

27ri 



X 



-2-1 



r(^ + i)nLi(^-«fc) 

F(z + n + 2 + /3) 



dz 



1 
27ri 



-2-1 



r(^ + i)nLi(^-«fe) 

F(z + n + 2 + /3) 



Rj, 



(268) 

we have that as — )• 00, 



where S-^ is the same as in (254). If x G [0, 1], by Lemma I h) 

the contour integral over vanishes for any x 7^ 0, and from ( |25"4 ), we find that "^JLi 

converges for any x G C. Thus the contour integral over S in (239) is a power series in x for all 
X G (0,1). 

To show that the right-hand side of (244) is a polynomial in x of degree n — 1, we apply the 
identity that for < x < 1 and k £ Z, 



c-^-iF(z + A; + l) 
27ri /s F(z + A; + 2 + /3) 



r(i + /3) 



(269) 
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The proof of (269) is similar to that of (255), Noting the similarity between the contour integral 
in (269) and that in (244) by arguments like above we have 



1 



2m Js T{z + k + 2 + P) 



dz 



j=o' 



-2-1 



res 



r{z + j + k + 2) 



Ef 

j=0 

oo T-rj-1 



-j-k-1 r{z + k + 2 + (3) UiZii^ + k + l) 



i-j + 2 + m 



(270) 



^ r(l + /3)j! ^ 



j=0 



Comparing (270) with the expansion 



(271) 



we prove (269 ). 



Similar to (257), there are coefficients cg , ci , . . . , Cn such that 

r(^+i)nLi(^-«^-) nz+k+i) 



T{z + n + 2 + P) 



k=0 



Viz + k + 2 + 13)'' 



where Cn 



In (272), the value of c„ is obtained, when /3 Z, by substituting z 
(272) becomes 

n 

r(-n - /3) J](-n - 1 - /3 - afc) = c„r(-/3). 



nLi(^+/3) • 

(272) 

— (n + 1 + /?), so that 



k=l 



and further by THopital's rule when /3 G Z. Summarizing the results obtained, we have 



27ri 



-13 



-2-1 



r(^ + i)nLi(^-«fe) 

r(z + n + 2 + /3) 



dz 



x\ 



k=0 



(273) 



(274) 



and hence prove that the right-hand side of (244) is a monic polynomial of degree n. 
We need to show that the orthogonal relation ( |241 ) holds. For m = 1, . . . , n. 



X 



-2-1 



dz x'^'^dx 



lim 

e-S>0+ 



1-e 



1-e 



r(^ + i)nLi(^-«fc) 

r(z + n + 2 + /3) 

X— ir(z+i)nLi(^-«fc: 

r{z + n + 2 + 13) 

r(^ + i)nLi(^-«fc) 



X ^x°'"'dx 



r{z + n + 2 + f3) 

+ 1) nfc=l,...,n(2 - "fc) 

^^!!! dz 

2 r(z + n + 2 + /3) 



1 x°''^dx 
dz 



(275) 



From the formula (244), we see that if we can show the last contour integral in (275) vanishes, 
then the proof is done. To see that, we take a = 6 = C, a large positive number, in the definition 
of the contour S in (250). By Lemma I \h) the integrand in the last contour integral in (275) 
is 0{\z\~^'^^^^). Thus the contour integral itself should be 0{C~^^~^^'^). Let C — )• oo, we prove 
that the contour integral vanishes. 
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B Construction of contours 

First we construct S of Figure [5] Let S be of the shape (see Figure [s]) 



S := Si U U Sa, where < 
Note that for a large enough a, we have the fohowing 



Si 
S3 



xo + (1 — + t, t < —a, 
xq + it, —a <t<a, 

Xo + (1 + — t, t > a. 



(276) 



Conditions for a 

1. For 2; G Si U S3, RF{z) decreases as z moves to the left. 

2. For t G {a, 00) U (-00, -a), ^F{it) < 3?F(xo). 



We assume the parameter a in (276) satisfy the two conditions above. 



To show that ^F{z) attains its maximum on S at xq, it then suffices to show that 'RF{z) 
decreases as z moves along S2 upward (resp. downward) from xq to xq ± ia. We need only to 
consider the case z G S2 H C+ due to the symmetry of KF(z) about the real axis. To this end, 
since -F'(xo) = F"{xq) = 0, it suffices to show that for all i > 0, 



d'^{^F{xQ + it)) 



-KF"(xo + it) 



1 



1 



1 



Xo + it (xo — a + ity xq — b + itY 
XQ , (xq - af - t^ , (xq - hf - 



+ 



+ 



t' 



(277) 



x2 + t2 ((xo-a)2 + t2)2 5)2+^2)2 



< 0. 



Using the expressions of xq, xo — a and xq — 6 in (172) and (173) and taking the substitution 

* = 2(l+r3) ^' ^^'^^ 

An _|_ ^3^2 

;(s), (278) 

yi -r r-)- 

where 



2n3 



G{s) 



(1+r^) 



(l+r2 
1 



1 — r s 



2o2 



(279) 



(l + r2)2+4(l+r3)2s2 (l + s2)2 (l + r2>,2)2- 

We need to show G[s) > when s > 0. Since G(0) = 0, it suffices to show G'{s) > for s > 0. 
Writing 



G'{s) 



[(l + r2)2 + 4(l + r3)2s2](i + s2)2 



X [(3(1 + r2)2 - 4(1 + r^f) + ((1 + r'^f + 4(1 + r'^f)s^] . (280) 

We check that the coefficient 

3(1 + r'^f - 4(1 + r^f = {-2r^ + VZr"^ + \/3 - 2)(\/3(l + r^) + 2(1 + r^)) (281) 

is positive for r G (— l,ro), hence G'{s) > for s > given that r G (— l,ro). Therefore our 
construction of S is valid. 

Next we want to construct F^ and F^ of Figure [5] such that 3?-F(z) attains its minimum on 
Fb (resp. Fa) near xq. Since F{^) = F{z), we consider only the parts of Ff, and F^ on the upper 
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half plane, and let the lower parts of Th and Ta be obtained by reflection. We construct first 



^prelim 



(resp. pP'''^^™) as an approximation of F), n C+ (resp. Ta H C+), and then deform 



prelim 



(resp. pP'''^^™) locally around xq to obtain n C+ (resp. Ta H C+). 

In our construction of the contours, we need a technical result that F'{z) has no zeros in 
the first quadrant (except for xq on the boundary). To see it, we consider the contour integral 
^ !f(j ^ F"{z)/F'{z)dz, where Cquad is a quarter-circle contour in the first quadrant, centered 



at and locally deformed at a, b and xq to exclude the singularities, see Figure [9j The value of 
this contour integral is the number of zeros of F'{z) enclosed by Cquad- By direct calculation, 
we see that the contour integral vanishes as the radius of Cquad approaches oo, which implies 
that F'(z) has no zero in the first quadrant. 




Xq + ia 



ta 




a Xq 



Figure 8: The actual shape of E. Figure 9: The contour Cquad- 

First we construct pP''°^™_ By the local property of F{z) around xq, we have that the 
direction f is a steepest-ascent direction of ?fiF{z) at xq. we construct the gradient flow line 
of the vector field S/^F(z) on the complex plane, (with the complex plane identified with M^,) 
starting from xq with initial direction j. The value of K-F(z) increases as z moves along the 
flow line. Consider the region 

B = {z €C+\0 <Qz < Ni, < "Siz < N2 and z is to the right of S,} (282) 

Here we choose A^i and such that they satisfy the following 



Conditions for A^i and A''2 

1. As z moves along the horizontal line M -|- iNi to the right, "SiFlz) is increasing. 

2. For z in the vertical line between N2 and + ^-^ii ^F{z) > ^F{xo). 

It is not difficult to check that if A'^i is large enough, then Condition [T] is satisfied, and for any 
fixed A^i, if N2 is large enough, then Condition [2] is satisfied. So the region B is well defined. 
We denote by the part of the gradient flow line in B and connected to xq. Since along the 
gradient flow line ?fiF{z) is eventually unbounded, F^ has to hit the boundary of B. The value 
of ?R-F{z) on S is less than ^F{xo), so F^ cannot hit S. For the same reason it cannot hit 
the line segment between xq and b, or the imaginary axis. Then there are only two cases to 
consider: 



1. Tb hits the ray (6, -|-oo): We define r^'^'^'™ as F^, as shown in Figure 



10 
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2. Tb hits the top or the right side of B: We define pP'''^'™ as the combination of and 
the right side (and possibly also the top) of B between its intersection with Tb and its 
intersection with the real axis, as shown in Figure 11 




Figure 10: One possible construction of Figure 11: Another possible construction 

-ppreiim ^j-^g^^ ]^j^g ^]-^g bottom of i?. of pP''*'^™ when F^ hits the right side of 

B, and then pP''*'^™ the union of F^ and 
the part of the right side of B (in boldface) 
between F^ and the real axis. 

It is clear that in either case, 3?-F(z) has xq as the unique minimum on pP'''^'™, 
The construction of pP'''^''™ ig done in a similar way. Denote 

^ = {z G C+ I z is to the left of S and < 9z < a}, (283) 

and construct the gradient flow line of V5R-F(z) starting from xq with initial steepest ascent 
direction We denote the part of the flow line in A and connected to xq by F^i. By the same 
argument as for F^, we have the following two possible cases: 
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1. F^ hits (0, a): we define Pa'^'^^™ as F^, as shown in Figure 

2. F^ hits the left side of A, i.e., the vertical line between and ia: we define Pa'^'^^'™ as the 
union of F^ and the vertical line between and the intersection of dA and F^, as shown 



in Figure 13 




Figure 12: One possible construction of Figure 13: Another possible construction 

■ppreimi ^j^gj^ p^ j^j^g (Q^a). of pg'''''™ when F^ hits the left side of B, 

and then pP''*'^™ the union of F^ and the 
vertical line (in boldface) between F^ and 
the real axis. 

Note that in case [l| ^F{z) has xq as the unique minimum on pP''*'^™ due to the property of 
the flow line F^, but in case [2], this is not always true. If a is big enough, or equivalently the 
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parameter r £ (— l,ro) defined in (171) is away from rg, we have that for z = it (t > 0) on the 
positive imaginary axis 



^ ^ ^ ^ = -QF'iit) = r + TT T--< 0, (284) 

t ' a t ' b 

and then K-F(z) is decreasing as z = it moves from upward along the imaginary axis. Hence 
the definition of pP^*^^™ ig valid if (284) holds for all t > 0. Numerical result shows that (284) 



holds if a > 0.3436, or equivalently r < —0.7150. 

Remark 11. Although more sophisticated construction of ra"^*^^™ allow smaller a, there is an 
essential difficulty if we let a — )• 0"''. Since we require Ta (and so pP'''^^™) not to intersect with 
(—00,0), there is zq S r^'^*'^™ n (0, a). We need ^F{zo) > KF(xo). However, it is not difficult 
to show that if a is very close to 0, or equivalently, r is very close to ro, 

^F{x) < RF{xo), for all x e (0, a). (285) 

and then the construction of pP'''^''™ ig impossible. Numerical result shows that ( 285 ) holds if 



a < 0.1184, or equivalently r > —0.6827. To consider the case that a is very close to 0, we need 
to allow Ta to cross (— oo,0), and radically change the contour integral formula. We do not 
pursue it in this paper. 

The contours pP''^^™ and pP'''^'™ with their reflections in the lower half plane are still not 
satisfactory contours, because (1) they intersect with S at xq, (2) they are not described by 
explicit formulas around xq, the point around which we do steepest-descent analysis. To solve 
these defects, we first choose a small enough e > and deform the part of p^^^^^'™ (resp. pP''''''™) 

that is within distance e to xq into the straight line segment between xq to xq + e~e (resp. 
between xq to XQ + e^e) such that "SiFlz) still increases as z moves along the contours away from 
Xq. Then we deform the line segment between xq and xq + e~ c^^n~* (resp. the line segment 
between xq and xq + c^^n~i) into the vertical line segment between xq — 

Xq + e 4 Cg n 4 (resp. the vertical line segment between xq + '^(^z * ^-^id xq + e 4 Cg n 4 . 
Then we obtain Ta n C+ (resp. F;, n C+) and the final contour Ta (resp. F;, is obtained by 
reflection about the real axis. 
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